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This publication describes and illustrates the use 
of the IBM FORTRAN IV language for the IBM Time Sharing 
System (TSS>. 

The IBM FORTRAN IV language is a symbolic 
programming language that parallels the symbolism and 
format of mathematical notation. It provides 
programming features and facilities that can be used in 
FORTRAN programs to solve mathematical problems. 

The material in IBM FORTRAN IV is arranged to 
provide a quick definition and syntactical reference to 
the various elements of FORTRAN IV by means of a box 
format. Each element is described, with appropriate 
examples of possible use. 

The reader should have a basic knowledge of the 
FORTRAN language. While some information relating 
FORTRAN IV to TSS is presented, most of the necessary 
guidance required by a FORTRAN user to perform a task 
is given in IBM Time Sharing System: FORTRAN 
Programmer's Guide , GC28-2025. 



PREFACE 



This manual describes the IBM Time 
Sharing System FORTRAN IV language. The 
material in this publication is arranged to 
provide a quick definition and syntactical 
reference to the various elements of 
FORTRAN IV by means of a box format. Each 
element is described, with appropriate 
examples of possible use. 

Five appendixes give additional 
information for writing a FORTRAN IV source 
program : 

A: FORTRAN Comparison 

B: Table of Source Program Characters 

C: Other FORTRAN Statements Accepted by 

ISM FORTRAN IV 
D: FORTRAN- Supplied Subprograms 
E: Examples of FORTRAN- Written Programs 



The user should have a basic knowledge 
of the FORTRAN language before using this 
publication. While some information 
relating FORTRAN IV to TSS is presented , 
most of the necessary guidance required by 
a FORTRAN user to perform a task is given 
in IBM Time Sharing System; FORTRAN 
Programmer's Guide , GC28-2025. 
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INTRODUCTION 



IBM TIME SHARING SYSTEM FORTRAN IV 

TSS FORTRAN IV includes a language, a compiler, and a set of system- 
supplied subprograms. 

For more information regarding the compiler and system- supplied pro- 
grams, see FORTRAN Programmer's Guide and FORTRAN IV Library Subpro- 
grams . The FORTRAN compiler operates under the control of TSS, which 
provides the compiler with input/output and other services. Object pro- 
grams generated by the compiler operate under TSS control and are depen- 
dent on it for similar services. 

The TSS FORTRAN IV language is compatible with and encompasses the 
American National Standards Institute (ANSI) FORTRAN, including its 
mathematical subroutine provisions. Source programs written in FORTRAN 
consist of a set of statements constructed from the elements of the lan- 
guage described in this publication. 

SPECIAL FEATURES OF TSS FORTRAN IV 

TSS FORTRAN IV is a further development of previously implemented 
FORTRAN systems and contains many of the features of these systems (See 
"Appendix A: FORTRAN Comparison"). Ttie following features facilitate 
the writing of source programs and reduce the possibility of coding 
errors: 

1. Variable Attribute Control : Ihe attributes of variables and arrays 
may now be explicitly specified in the source program. This faci- 
lity is provided by a single explicit specification statement that 
allows a programmer to: 

a. Explicitly type a variable as integer, real, complex, double 
precision, or logical. 

b. Specify the number of storage location bytes to be occupied by 
each variable or member of an array. 

c. Specify the dimension of an array. 

d. Specify data initialization values for variables. 

2. Adjustable Array Dimensions : The dimensions of an array in a sub- 
program may be specified as variables; when the subprogram is 
called, the absolute array dimensions are substituted. 

3. Additional Format Code : An additional format code - G - can be 
used to specify the format of numeric and logical data. Previously 
implemented format codes are also permitted. 

4. Mixed Mode : Expressions may consist of constants and variables 
that are of the same and/or different types and lengths. 

5 . Namelist I/O : Formatting of input/output data is facilitated by 
reading and writing operations without reference to a FORMAT state- 
ment or list. 

6. Spacing Format Code : Ihe T format code allows input /output data to 
be transferred beginning at any specified position. 

7. Literal Format Code : Apostrophes may be used to enclose literal 
data. 
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ELEMENTS OF THE LANGUAGE 



STATEMENTS 

FORTRAN statements are composed of FORTRAN key words used in conjunc- 
tion with the basic elements of the language (constants, variables, and 
expressions). The five categories of FORTRAN statements are: 

1. A rithmetic and Logical Assignment Statements ; Replace the current 
value of a designated variable after calculations have been 

performed. 

2 . Control Statements : Govern the flow and terminate the execution of 
the object program. 

3. Input /Out put Statements : Exchange information between a user's 
program and a named collection of data, called a data set. 

4 . Specification Statements : Declare the properties of variables, 
arrays, and subprograms (e.g., type and amount of storage reserved) 
and describe the format of input or output data, 

5. Subprogram Statements : Define and name functions and subroutines. 

CODING FORTRAN STATEMENTS 

Card Input 

The statements of a FORTRAN source program can be written on a stan- 
dard FORTRAN coding form. Form X28-7327 (see Figure 1). FORTRAN state- 
ments are written one to a line from columns 7 through 72. If a state- 
ment is too long for one line, it may be continued on as many as 19 suc- 
cessive lines by placing any character, other than a blank or zero, in 
column 6 of each continuation line. For the first line of a statement, 
column 6 must be blank or zero. 

Columns 1 through 5 of the first line of a statement may contain a 
statement number consisting of from one through five decimal digits. 
Leading zeros in a statement number are ignored. Statement numbers may 
appear anywhere in columns 1 through 5 and may be assigned in any order; 
the value of statement numbers does not affect the order in which the 
statements are executed in a FORTRAN program. 

Columns 73 through 80 are not significant to the FORTRAN compiler and 
may, therefore, be used for program identification, sequencing, or any 

other purpose. 

Comments to explain the program may be written in columns 2 through 
80, if the letter C is placed in column 1. Comments may appear anywhere 
within the source program except immediately preceding a continuation 
line. They are not processed by the FORTRAN compiler, but are printed 
on the source program listing. Blanks may be inserted where desired to 
improve readability. 

Keyboard Input 

The conversational capability of TSS FORTRAN allows statements to be 
entered via a keyboard at a terminal. The rigid card column rules are 
relaxed for this means of input. A detailed description of keyboard 
input is contained in FORTRAN Programmer's Guide . 




Figure 1. FORTRAN coding form 



CONSTANTS 

A constant is a fixed, unvarying quantity. The three classes of con- 
stants are those that deal with numbers (numerical constants), truth 
values (logical constants), and literal data (literal constants) . 

Numerical constants may be .integer, real, or complex numbers? logical 
constants may be .TRUE. or . FALSE. ; literal constants may be a string 
of alphameric and/or special characters enclosed by apostrophes. 



INTEGER CONSTANTS 



Definition 



h~ 






| Integer Constant - a whole number written without a decimal point; 
| it occupies two or four storage locations; four storage locations is 
| the standard length. 



I 



Maximum Magnitude: 2147483647, i.e., (2 3 *~1) 



An integer constant may be positive, zero, or negative; if unsigned, 
it is assumed to be positive. Its magnitude must not be greater than 
the maximum and may not contain embedded commas. 

Examples ; Valid Integer Constants: 



91 

173 

-2147483647 

-12 
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Invalid Integer Constants: 



0.0 

27. 

3145903612 

5,396 



(contains a decimal point} 
(contains a decimal point) 
(exceeds the allowable range} 
(contains embedded comma} 



REAL CONSTANTS 



Definition 



y~ 



4 



R eal Constant — a number with a decimal point optionally followed 
by a decimal exponent, or an integer constant followed by a decimal 
exponent. The exponent may be wri' .^n as the letter E or D followed 
by a signed or unsigned, one- or two-digit integer constant. A real 
constant may assume one of two forms: 

1. From 1 through 7 decimal digits, optionally followed by an E and 
a decimal exponent. This form occupies 4 storage locations. 

2. Either 1 through 7 decimal digits, followed by a D and a decimal 
exponent, or 8 to 16 decimal digits, optionally followed by a D 
and a decimal exponent. This form occupies 8 storage locations 
and is sometimes referred to as a double- precis ion constant. 



Magnitude: 






(either form) 
ately 10"). 



or 16-" through 16° 3 (i.e., approxim- 



— j 



A real constant may be positive, zero, or negative (if unsigned, it 
is assumed to be positive) and must be of the allowable magnitude. It 
may not contain embedded commas. Hie decimal exponent E or D permits 
the expression of a real constant as the product of a real constant 
times 10, raised to a desired power. If a decimal exponent is given, 
the decimal point is not required. 

Examples : Valid Real Constants (4 storage locations): 

+ 0. 

-999.9999 

0.0 

5764.1 

7.0E+0 

19761. 25E+1 

7E3 

7.E3 

7.0E3 

7.0E03 

7.0E+0 3 

7.0E-0 3 



(i.e. , 7.0 x 10° = 7. 0) 

(i.e., 19761.25 x 10* = 197612.5) 



(i.e., 7.0 x 103 = 7000.0) 



(i.e. , 7. x 10- 3 = . 007) 



Valid Real Constants (8 storage locations): 



21.987538 29457168 

1.0000000 

7.9D3 

7.9D03 

7.9EH-03 

7.9EH3 

7.9D-03 

7.9D0 

0.0D0 

7D3 



(i.e., 7. 9 x 10 3 = 7900.0) 



(i.e., 7.9 x 10- 3 = .0079) 

(i.e. , 7.9 x 10° = 7. 9) 

(i.e. , 0. x 10° = 0. 0) 

(i.e. , 7 x 10 3 = 7000) 



Invalid Real Constants: 

(missing a decimal point) 

3 f H71.1 (contains embedded comma) 

l.E (missing a one- or two-digit integer 

constant following the E; note that it is not 

interpreted as 1,0 x 10°) 
7 .9D (missing a 1- or 2-digit integer 

constant following the D) 
1.2E+113 (F is followed by a 3-digit 

integer constant) 
21.3D90 (value exceeds the magnitude permitted; 

that is, 21.3 x 10»°>16« 3 ) 
23.5E+97 (value exceeds the magnitude permitted; 

that is, 23.5 x 10* 7 >16 63 ) 



COMPLEX CONSTANTS 

r — ' * ■ * ' ' "' — ■— ~~ 

| Definition 

, .^ . _ ._ __ _„_ _„___„ 

Complex Constan t — an ordered pair of signed or unsigned real con- 
stants separated by a comma and enclosed in parentheses. These real 
constants may assume one of two forms: 



1. From 1 through 7 decimal digits, optionally followed by an F 
decimal exponent. In this form, each number in the pair occu- 
pies 4 storage locations. 

2. Either 1 through 7 decimal digits, followed by a D decimal 
exponent, or 8 through 16 decimal digits, optionally followed by 
a D decimal exponent. In this form each number in the pair 
occupies 8 storage locations. 

Magnitude: (either form) or 16-" through 16* 33 (i.e., approxim- 
ately 10 7S ) for each real constant in the pair. 



The real constants in a complex constant may be positive, zero, or 
negative (if unsigned they are assumed to be positive), but they must be 
within the given range. The first real constant in a complex constant 
represents the real part of the complex number; the second represents 
the imaginary part of the complex number. 

Examples : Valid Complex Constants: 

(3.2,- 1.86) (has value 3.2-1.861) 

(-5.0E+03,.16E4-02) (has value - 5000. +16. 0i> 

(4.0E+03,.16E*02) (has value 4000. +16. Oi) 

(2.1,0.0) (has value 2.1+0. 01) 

(4. 1D+2,1. 9736148) (has value 470. +1. 97361481) 

where i 

Invalid Complex Constants: 

(292704,1.697) (real part does not 

contain decimal point) 
(1.2E113,279.3) (real part contains 

invalid decimal exponent) 
C .0034E4, .005D6) (parts differ in length) 
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LOGICAL CONSTANTS 



j Definition j 

^ f 

| Logical Constant — The two logical values are: | 

I ! 

| .TRUE, I 

I -FALSE- | 



A logical constant mist be preceded and followed by a period. The 
logical constants .TRUE. and .FALSE. specify that the value of the 
logical variable they replace , or the term of the expression they are 
associated with, is true or false, respectively (see "Logical 

Express ions"). 

LITERAL CONSTANTS 



•1 



| Definition 

Y f 

| Literal Constant — a string of alphameric and/or special characters | 
| enclosed in apostrophes. j 

L J 

Ihe minuter of characters in the string, including blanks, may not be 
greater than 255. Since apostrophes delimit literal data, a single apo- 
strophe within such data is represented by double apostrophes. An 
alternative form for a literal constant is wH immediately followed by a 
string, of length w, of alphameric and/or special characters. A single 
apostrophe within such data is represented as a single apostrophe. 

Examples ; 

•DATA 1 

•INPUT /OUTPUT AREA NO. 2* 

• X-COORDI NATE Y-COORDINATE Z-COORDINATE" 

•3. 14' 
•DON' *T 9 
5HDON*T 



HEXADECIMAL CONSTANTS 

f ^ 

I Definition | 

j. 4 

[ Hexadecimal Constant — the character Z followed by a number formed | 
| from the set through 9 and A through F. | 

L . . . J 

Hexadecimal constants may be used only as data initialization values. 

One storage location contains two hexadecimal digits. If a constant 

is specified as an odd number of digits, a leading hexadecimal zero is 

added to fill the storage location. The internal form of each hexade- 
cimal digit is: 

0-0000 **-0100 8-1000 C-1100 

1-0001 5-0101 9-1001 D-1101 

2-0010 6-0110 A-1010 E-1110 

3-0011 7-0111 B-1011 F-llll 



Examples : 

Z1C49A2F1 
Z BADFAD 

The maximum number of digits allowed in a hexadecimal constant 
depends upon the length specification of the variable being initialized 
Csee "Variable Types and Length Specifications"). The following list 
shows the maximum number of digits for each length specification. 

Length Specification of Variable Maximum Number of Hexadecimal Digits 
16 32 

8 16 

4 8 

2 t* 

1 2 

If the number of hexadecimal digits is greater than the maximum, the 
leftmost digits are truncated; if the number of digits is less than the 
maximum, leading hexadecimal zeros are added. 



SYMBOLIC NAMES 



j Definition j 

v , 

J Symbolic Name — from 1 through 6 alphameric (i.e., through 9, or 1 
| alphabetic, A through Z and $) characters, the first of which must j 
j be alphabetic. | 



Symbolic names are used in a program unit (i.e., a main program or a 
subprogram) to identify elements in the following classes. 

• An array and the elements of that array tsee "Arrays") 

• A variable (see "Variables") 

• A statement function (see "Statement Functions") 

• An intrinsic function (see Appendix D) 

• A FUNCTION subprogram (see "FUNCTION Subprograms") 

• A SUBROUTINE subprogram (see "SUBROUTINE Subprograms") 

• A block name (see "BLOCK DATA Subprograms") 

• An external procedure that cannot be specified as either a SUBROU- 
TINE or FUNCTION subprogram (see "The EXTERNAL Statement" ) 

Symbolic names must b€> unique within a class in a program unit and 
can identify elements of only one class with the following exceptions. 

A block name can also be an array, variable, or statement function 
name in a program unit. 

A FUNCTION subprogram name must also be a variable name in the FUNC- 
TION^ subprogram. 

Once a symbolic name is used as a FUNCTION subprogram name, a SUBROU- 
TINE subprogram name, a block name, or an external procedure name in any 
unit of an executable program, no other program unit of that executable 
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The three ways a programmer may declare the type of a variable are by 
use of the: 

1. Predefined specification contained in the FORTRAN language, 

2. IMPLICIT specification statement, 

3. Explicit specification statements. 

The optional- length specification of a variable may be declared only 
by the IMPLICIT or Explicit specification statements. If, in these sta- 
tements, no length specification is stated, the standard length is 
assumed (see "The Type Statements •) . 



TYPE DECLARATION BY THE PREDEFINED SPECIFICATION 

The predefined specification is a convention used to specify 
variables as integer or real: 

1. If the first character of the variable name is I, J, K, L, M, or N, 
the variable is integer of standard length, **. 

2. If the first character of the variable name is any other letter, 
the variable is real of standard length, **. 

This convention is the traditional FORTRAN method of implicitly spe- 
cifying the type of a variable as being either integer or real. In all 
examples that follow, it is presumed that this specification holds, 
unless otherwise noted. 



TYPE DECLARATION BY THE IMPLICIT SPECIFICATION STATEMENT 

This statement allows a programmer to specify the type of variables 
in much the same way as was specified by the predefined convention. 
That is, in both, the type is determined by the first character of the 
variable name. However, the programmer, using the IMPLICIT statement, 
has the option of specifying which initial letters designate a particu- 
lar variable type. Further, the IMPLICIT statement is applicable to all 
types of variables — integer, real, complex, and logical. 

The IMPLICIT statement overrides the variable type as determined by 
the predefined convention. For example, if the IMPLICIT statement spe- 
cifies that variables beginning with the letters A through M are real 
variables, and variables beginning with the letters N through Y are 
integer variables, then the variable ITEM (which would be treated as an 
integer variable ander the predefined convention) is now treated as a 
real variable. Note that variables beginning with the letters Z and $ 
are (by the predefined convention) treated as real variables. The 
IMPLICIT, statement is presented in greater detail in the section "Type 
Statements ." 

TYPE DECLARATION BY EXPLICIT SPECIFICATION STATEMENTS 

Explicit specification statements (INTEGER, REAL, COMPLEX, and LOGIC- 
AL) differ from the first two ways of specifying the type of a variable, 
in that an explicit specification statement declares the type of a par- 
ticular variable by its name, rather than as a group of variables begin- 
ning with a particular character. 

For example, assume: 
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1. That an IMPLICIT .specif icati on statement overrode the predefined 
convention for variables beginning with the letter I # by declaring 
them to be real. 

2. That a subsequent Explicit specification statement declared that 
the variable named ITEM is complex. 

Then, the variable ITEM is complex and all other variables beginning 
with the character I are real. Note that variables beginning with the 
letters J through N are specified as integer by the predefined 
convention. 

The Explicit specification statements are discussed in greater detail 
in the section "Type Statements. " 

ARRAYS 

A FORTRAN array is a set of variables identified by a single variable 
name. A particular variable in the array may be referred to by its 
position in the array (e.g., first variable, third variable, seventh 
variable, etc.). Consider the array named NEXT, which consists of five 
variables, each currently representing these values: 

273, 41, 8976, 59, and 2 



NEXTC1) 
NEXTC2) 
NEXT (3) 
NEXT(*4) 

NEXT (5) 



represents 2 73 
represents Ml 
represents 8976 
represents 59 
represents 2 



Each variable in this array consists of the name of the array (i.e., 
NEXT) followed by a number enclosed in parentheses, called a subscript. 
The variables that constitute the array are called subscripted 
variables. Therefore, the value of the subscripted variable NEXT(l) is 
273; the value of NEXT(2), is i»lj etc. 

The subscripted variable NEXT(I) refers to the " Ith " subscripted 
variable in the array, where I is an integer variable that can assume a 
value of 1, 2, 3, 4, or 5, in this array. 

To refer to the first element of an array, the array name must be 
subscripted; the array name does not represent the first element. The 
number of subscripts must correspond to the declared dimensionality, 
except in the EQUIVALENCE statement, which is explained in "EQUIVALENCE 
Statement • " 

SUBSCRIPTS 

A subscript Is an index that specifies one of the coordinates that 
identify a particular element of an array. From one to seven subscripts 
are used in accordance with the dimensionality of the array being sub- 
scripted. Multidimensional subscripts are separated by commas. The 
subscripts, enclosed in parentheses, follow the array name. 
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I General Form 

| . 

Subscripts — may be one of seven forms 2 

v 

c' 

v*c f 

v-c f 

c*"v 

c*v+c' 

c*v-c • 

where v represents an unsigned, nonsubscripted f integer variable. 

c and c' represent unsigned integer constants. 






Whatever subscript form is used, its evaluated result must always be 
greater than zero. For example, when reference is made to the sub- 
scripted variable VtI-2), the value of I should be greater than 2. 

Examples 1 

ARRAY (IHOLD) 
NEXT (19) 
MATRIX (1-5) 
A(5*L) 
W(4*M+3) 

These are valid subscripted variables for their corresponding arrays: 



Array Name 

A 

TABLE 

B 

MATRIX 



Subscripted Variable 

A(5, 100, J, K+2) 

TABLE (1, 1, 1, 1, 1, 1, 1) 

B(I, J, K, L, M, N) 

MATRIX( 1+2 , 6*JOB- 3,KFRAN) 



Consider the following array, named LIST, consisting of two subscript 
parameters, the first ranging from 1 through 5, the second from 1 
through 3, 



Column 1 


Column 


2 


Column 3 


82 


4 




7 


12 


13 




14 


91 


1 




31 


2** 


16 




10 


2 


8 




2 



Row 1 

Row 2 

Row 3 

Row 4 

Row 5 

Hie correct reference for the number in row 2, column 3 f would be 

LIST (2,3) 

LIST (2,3) has the value 14; LIST (4,1) has the value 24. 

Ordinary mathematical notations might use LIST i,j to represent any 
element of the array LIST. In FORTRAN, this is written as LIST (1 , J) t 
where I equals 1, 2, 3, 4, or 5, and J equals 1, 2, or 3. 

As a further example, consider the array named COST, consisting of 
four subscript parameters. Ttiis array might be used to store all the 
premiums for a life insurance applicant, given (1) age, (2) sex, (3) 
health, and (4) size of life insurance coverage desired. A code number 
could be developed for each statistic; IAGE represents age; ISEX, sex; 
IHLTH, health; and ISIZE, policy size desired (see Table 1). 
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Table 1. Insurance premium codes 

r , ._. 



AGE 



k' 



Age in years 

1-5 

6-10 



96-100 



Code 

IAGE=1 
IAGE=2 



IAGE=2 






HEALTH 



Health 

Poor 

Fair 
Good 

Excellent 



Code 

JHLTH=1 
IHLTH=2 
IHLTH=3 
IHLTH=4 



\ 

I 

V 

I 



SEX 



4 



Sex 




Code 




Male 




ISEX=1 




Femal 


e 


ISEX=2 


4 








1 




POLICY 


SIZE 


— \ 






Dollars 




Code 




1,000 




ISIZE=1 




3,000 




ISIZE=3 




5,000 




isizE=a 




10,000 




ISIZE=5 




25,000 




ISIZE=6 




50,000 




ISIZE=7 




100,000 




ISIZE=8 








Suppose an applicant is 14 years old, male, in good health, and 
desires a policy of $25,000, From Table 1, the statistics can be repre- 
sented by these codes 



IAGE=3 
ISEX=1 

1HLTH=3 
ISIZE = b 



Cll - 15 years old) 

(male) 

(good health) 

($25,000 policy) 



Thus, COST (3, 1, 3, 6) represents the premium for a policy, given 
the statistics above. Note that IAGE can vary from 1 to 20, ISEX from 1 
to 2, IHLTH from 1 to 4 , and ISIZE from 1 to 8. The number of sub- 
scripted variables in the array COST is the number of combinations that 
can be formed for different ages, sex, health, and policy size available 
- a total of 20x2x4x8 or 1280. Therefore, up to 1280 different premiums 
may be stored in the array named COST. 



DECLARING THE SIZE OF AN ARRAY 

The size of an array is determined by the number of subscript parame- 
ters of the array and the maximum value of each subscript. This infor- 
mation must be given for all arrays before using them in a FORTRAN pro- 
gram, so that a storage area of sufficient size may be reserved. 
Declaration of this information is made by a DIMENSION statement, a COM- 
MON statement, or by one of the Explicit specification statements 
(INTEGER, REAL, COMPLEX, and LOGICAL); each is discussed in "Specifica- 
tion Statements." 

ARRANGEMENT OF ARRAYS IN STORAGE 

Arrays are stored in ascending storage locations, with the value of 
the first of their subscripts increasing most rapidly, and the value of 
the last increasing least rapidly. 

Examples : The array named A, consisting of one subscript parameter, 
^hich varies from 1 to 5, appears in storage as 

A(l) A(2) AC3) A(4) AC5) 
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The array named B consists of two subscript parameters? the first 
subscript varies over the range from 1 to 5, and the second varies from 
1 to 3. The array appears in ascending storage locations in this order 

B<1,1) B(2 f l) BC3.1) BU,1> B(5 f l> B(l,2) B(2,2) B<3,2) B(4, 2V— -| 



1 »*B<5,2) BC1/J) B<2,3) B(3,3) B<4,3> B(5,3) 

Note that B(l,2) and B(l,3> follow in storage B<5,1) and BC5 r 2) # 
respectively. 

The following list is the order of an array named C that consists of 
three subscript parameters; the first subscript varies from 1 to 3; the 
second, from 1 to 2; and the third , from 1 to 3. 

C(l,l,l) C<2 # 1,1> C(3,l,l> C(l,2,l> C(2 # 2,l) €(3, 2,1) C( 1,1,2 V— | 



I ^C(2 f l72) C13,1,2> C(l,2,2> C(2.2,2) C(3,2,2) C(l,l f 3) C<2,1. 3>— — | 

1 ^C(3 # l # 17T7l7273)~C(2 # 2 # 3) C(3,2,3) 

Note that C (1,1,2) and CCl f l f 3) follow in storage C(3 f 2 # l) and CC 3, 2, 2) , 
respectively. 



EXPRESSIONS 

Expressions in their simplest form consist of a single constant or 
variable. They may also designate a computation or show a relationship 
between two or more constants and/or variables. Expressions may appear 
in arithmetic and logical assignment statements and in certain control 
statements, 

FORTRAN provides two kinds of expressions; arithmetic and logical. 
The value of an arithmetic expression is always a number whose type is 
integer, real, or complex- However, the evaluation of a logical expres- 
sion always yields a truth value: • TRUE. or .FALSE.. 



ARITHMETIC EXPRESSIONS 

The simplest arithmetic expression consists of a single constant, 
variable, or subscripted variable of the type integer, real, or complex. 
If the constant, variable, or subscripted variable is of the type integ- 
er, the expression is in the integer mode. If it is of the type real, 
the expression is in the real mode, etc. 

Examples s 

E xpression Type of Quantity Mode of Expression 

"3 Integer Constant Integer of length 4 

I Integer Variable Integer of length 4 

Real Constant Real of length 4 

Real Variable Real of length 4 

j.j.-,^ Real Constant Real of length 8 

B(2*I) Real Variable Real of length U 

(2.0,5.7) Complex Constant Complex of length 8 

Complex Variable Complex of length 8 

(Specified as such in a 

Type statement) 



3.0 

A 

3.14D3 



C 
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In the expression B(2*I), the subscript (2*1), which must always 
represent an integer, does not affect the mode of the expression. That 
is # the mode of the expression is determined solely by the type of con- 
stant, variable, or subscripted variable appearing in that expression. 

More complicated arithmetic expressions containing two or more con- 
stants and/or variables may be formed by using arithmetic operators that 
express the computations to be performed. 

Arithmetic Operators 

The arithmetic operators are: 

Arithmetic Operator Definition 

** Exponentiation 

* Multiplication 
/ Division 

♦ Addi ti on 

Subtraction 

RULES FOR CONSTRUCTING ARITHMETIC EXPRESSIONS : These are the rules for 
constructing arithmetic expressions that contain arithmetic operators: 

1. All desired computations must be specified explicitly. That is, if 
more than one constant, variable, subscripted variable, or function 
reference (see " SUBPROGRAMS") appears in an arithmetic expression, 
they must be separated from one another by an arithmetic operator. 
For example, the two variables A and B will not be multiplied if 
written as 

AxB or AB or A«B 

If multiplication is desired, then the expression must be written 
as 

A*B or B*A 

2. No two arithmetic operators may appear in sequence in the same 
expression. For example, these expressions are invalid 

A*/B and A*-B 

However, in the expression, A*-B, if the - is meant to be a minus 
sign rather than the arithmetic operator designating subtraction, 
then the expression could be written as 

A*(-B) 

In effect, -B will be evaluated first, and then A will be multip- 
lied with it (for further uses of parentheses, see Rule 6). 

3. The mode of an arithmetic expression is determined by the type of 
the operands (where an operand is a variable, constant, function 
reference, or another expression) in the expression. Table 2 indi- 
cates how the mode of an expression that contains operands of dif- 
ferent types may be determined using the operators: +, -, *. /. 

Table 2 shows a hierarchy of type and length specification (see 
"Type Statements") that determines the mode of an expression. For 
example, complex data that has a length specification of 16, when 
combined with any other types of constants and variables, results 
in complex data of length 16. 



1*4 



Table 2. Determining the mode of an expression containing 
variables of different types and lengths 



♦ - * / 



, + + 



| IOT EGER 

| (2) 

| + 

| INTEGER 
| (4) 

I + 

(REAL 
I (4) 



| REAL 
| (8) 

I 

| COMPLEX 
| (8) 



(COMPLEX 
1 (16) 



INTEGER 
(2) 



Int eger 
(2) 



-+ 



Real 

my 

, + 



Real 
(8) 



•+- 



Complex 
(8) 



Complex 
(16) 



INTEGER j REAL 



(4) 



Integer 
(4) 



Integer | Integer 
(4) | (4) 

Real 
(4) 



Real 
(8) 






Complex 
(8) 






Complex 
(16) 



(4) 



Real 
(4) 



Real 
(4) 



Real 
(4) 



Real 
(8) 



Complex 
(8) 



Complex 
(16) 



REAL 
(8) 



Real 
(8) 



+ 



Real 
(8) 



Real 
(8) 



Real 
(8) 



-+ 



f 



Complex 
(16) 



Complex 
(16) 






COMPLEX 
(8) 



Complex 
(8) 



Complex 
(8) 



--+- 



Complex 
(8) 



Complex 
(16) 



Complex 
(8) 



Complex 
(16) 

L 



COMPLEX 
(16) 



Complex 
(16) 



Complex 
(16) 

Complex 

(16) 
4 

Complex 
(16) 



Complex 
(16) 

. 4 

Complex 
(16) 



Assume that the type of the following variables has been specified 
as 



Variable Names 
ROOT, E 
A f I f F 
C, D 



Type 

Real variables 
Integer variables 
Complex variables 



Length Specification 

4,8 

4,2,2 

16,8 



Then, the following examples illustrate how constants and variables 
of differing types may be combined using the arithmetic operators 
(♦. -, /r ♦> 



Expression 

ROOT* 5 

A+3 

C+2.9D10 

E/F+19 

C-18.7E05 

A/I-D 



Mode of Expression 
Real of length 4 
Integer of length 4 
Complex of length 16 
Real of length 8 
Complex of length 16 
Complex of length 8 



The arithmetic operator denoting exponentiation (i.e.,**) may only 
be used to combine the types of constants, variables, and sub- 
scripted variables shown in Table 3, 

Table 3. Valid combinations using the arithmetic operator ** 



| Base Exponent | 

(Integer or Real (either length)** Integer or Real (either length)) 



[Complex (either length) 



** Integer (either length) 



Assume that the types of the following variables are as specified, 
and that their length specification is standard. 



Variable Names 

ROOT,E 

A, I, F 

C 



Type 

Real variables 
Integer variables 
Complex variables 
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Then the following examples illustrate how constants and variables 
of differing types may be combined using the arithmetic operator 



Examples : 

Expression 

ROOT** (A* 2) 

C**A 

ROOT**I 

I**F 

7.98E21**ROOT 

ROOT**2.1E5 

A**E 



Type 

(Real**Integer ) 
(Complex* *Integer) 
(Real**Integer) 
(Integer**Integer > 
(Real**Real) 
(Real**Real) 
( Integer* *Real ) 



Result 

(Real) 

(Complex) 

(Real) 

(Integer) 

(Real) 

(Real) 

(Real) 



Order of Computation : Where parentheses are omitted, or where the 
entire arithmetic expression is enclosed within a single pair of 
parentheses, the order in which the operations are performed is as 
follows: 



Operation 

Evaluation of Functions (see 

" S u bpr ogr a m s " ) 
Exponentiation (**) 

Multiplication and Division (* and /) 
Addition and Subtraction (♦ and -) 



Hierarchy 

1st (highest) 

2nd 
3rd 
4th 



In addition, if two operators of the same hierarchy (with the 
exception of exponentiation) are used consecutively, the component 
operations of the expression are performed from left to right. 
Thus the arithmetic expression A/B*C is evaluated as if the result 
of the division of A by B was multiplied by C. 

For example, the expression 

(A*B/C**H-D) 



is evaluated in this order 

a. C**I Call the result X 

b. A*B Call the result Y 

c. Y/X Call the result Z 

d. Z + D Final operation 



( expon ent iat ion ) 
(mul tipl ica tion ) 
(division) 
(addition) 



For exponentiation the evaluation is from right to left, 
expression 

A**B**C 



Hi us, the 



is evaluated as 



a. B**C 

b. A**Z 



Call the result Z 
Final operation 



Use of Parentheses : Parentheses may be used in arithmetic expre- 
ssions, as in algebra, to specify the order in which the arithmetic 
operations are to be performed. Where parentheses are used, the 
expression within the parentheses is evaluated before the result is 
used. 

For example, the expression 

(B+((A*B) *C)+A**2) 

is evaluated in this order 
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a. (A+B) Call the result X 

b. (X*C) Call the result Y 

c. A**2 Call the result Z 

d. B+Y+Z Final operations 

7. Integer Division : When division is performed using two integers, 
the answer is truncated and an integer answer is given. For 
example, if 1=9 and J ==2, then the expression (I /J) would yield an 
integer answer of *4 after truncation. 

LOGICAL EXPRESSIONS 

The simplest form of logical expression consists of a single logical 
constant, logical variable, or logical subscripted variable, the value 
of which is always a truth value (i.e., either .TRUE. or .FALSE.). 

More complicated logical expressions may be formed by using logical 
and relational operators. These expressions may be in one of three 
forms . 

1. Relational operators combined with arithmetic expressions whose 
mode is integer or real. 

2. Logical operators combined with logical constants (.TRUE. and 
.FALSE.) , logical variables, or subscripted variables. 

3. Logical operators combined with either or both forms of the logical 
expressions described in items 1 and 2. 

Item 1 is discussed in the following section, "Relational Operators 1 *; 
items 2 and 3 are discussed in "Logical Operators." 

Relational Operators 

The six relational operators, each of which must be preceded and fol- 
lowed by a period, are: 

Relational Operator Definition 

.GT. Greater than (>) 

.GE. Greater than or equal to (>) 

.LT. Less than (<) 

.LE. Less than or equal to (<) 

.EQ. Equal to (=) 

.NE. Not equal to (*) 

The relational operators express an arithmetic condition which can be 
either true or false. Only arithmetic expressions whose mode is integer 
or real may be combined by relational operators. For example, assume 
that the types of these variables have been specified as 

Variable Names Type 
ROOT, E Real variables 

A, I, F Integer variables 

L Logical variable 

C Complex variable 

Then, the following examples illustrate valid and invalid logical 
expressions using the relational operators. 

Examples : Valid Logical Expressions Using Relational Operators! 

(RO0T*A).GT.E 

A.LT.I 

£**2.7.EQ. (5*ROOT4-*0 

57.9.LE.U.7+F) 

.5.GE..9*ROOT 

E.EQ.27.3D+0S 
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Invalid Logical Expressions Using Relational Operators: 

C.LT.ROOT (Complex quantities may never appear in logical 

expressions) 

C.GE. (2.7 r 5.9E3) (Complex quantities may never appear in logical 

expressions) 

L.EQ. (A+F) (Logical quantities may never be joined by rela- 

tional operators) 

E**2.EQ97.1E9 (Missing period immediately after relational 

operator) 

.GT.9 (Missing arithmetic expression before relational 

operator) 

Logical Operators 

The three logical operators, each of which must be preceded and fol- 
lowed by a period, are as follows, (A and B represent logical constants 
or variables, or expressions containing relational operators). 

Logical Operator Definition 

.NOT. .NOT. A - if A is .TRUE., then .NOT. A has the value 

.FALSE.; if A is .FALSE., then .NOT. A has the value 

.TRUE. 

.AND. A.AND.B - if A and B are both .TRUE., then A.AND.B 

has the value .TRUE. ; if either A or B or both are 
.FALSE., then A.AND.B has the value .FALSE. 

.OR. A.OR.B - if either A or B or both are .TRUE., then 

A.OR.B has the value .TRUE.? if both A and B are 
.FALSE. , then A.OR.B has the value .FALSE. 

Two logical operators may appear in sequence only if the second one 
is the logical operator .NOT.. 

Only those expressions which, when evaluated, have the value .TRUE, 
or .FALSE. may be combined with the logical operators to form logical 
expressions. For example, assume that the types of these variables are 
as specified. 

Variable Names Type 
ROOT, E Real variables 

A, I, F Integer variables 

L, W Logical variables 

C Complex variable 

Then, the following examples illustrate valid and invalid logical 
expressions using both logical and relational operators. 

Examples : Valid Logical Expressions: 

(ROOT*A.GT.A) .AND.W 
L.AND. .NOT.(I.GT.F) 
(E+5.9D2.GT.2+E) .OR.L 
. NOT. W. AND.. NOT. L 
L.AND. .NOT.W.OR.I.GT.F 
(A**F.GT.ROOT) .AND. -NOT. (I.EQ.E) 

Invalid Logical Expressions: 

A.AND.L (A is not a logical expression) 

.OR.W (.OR. must be preceded by a logical 

expression) 
NOT.(A.GT.F) (missing period before the logical operator 

.NOT. ) 
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(C.EQ.I).AND.L 

L.AND..OR.W 
.AND.L 



(a complex variable may never appear in a log- 
ical expression) 

(the logical operators .AND, and .OR. must 
always be separated by a logical expression) 
(.AND. must be preceded by a logical 
expression) 



Order of Computations in Logical Expressions : Where parentheses are 
omitted, or where the entire logical expression is enclosed within a 
single pair of parentheses, this is the order in which the operations 
are performed. 



Operation 

Evaluation of Functions 

Exponentiation (*♦) 

Multiplication and division (* and /) 

Addition and subtraction (♦ and -) 

• LT • , «LE. , • EQ • , . NE. , • GT. , • GE. 

.NOT, 

.AND. 

.OR. 

For example, the expression 

(A.GT.D**B.AND..NOT.L.OR.N) 

is evaluated in this order 



Hierarchy 

1st (highest) 

2nd 

3rd 

4th 

5th 

6th 

7th 

8th 



1. 


D**B 


Call the result 


W 


2. 


A.GT.W 


Call the result 


X 


3. 


.NOT.L 


Call the result 


y 


a. 


X.AND.Y 


Call the result 


z 


5. 


Z.OR.N 


Final operation 





( exponent ia tion ) 
(relational operator) 
(highest logical operator) 
(second highest logical operator) 



Use of Parentheses in Logical Expressions ! Parentheses may be used in 
logical expressions to specify the order in which the operations are to 
be performed. Where parentheses are used, the expression contained 
within the innermost pair of parentheses is evaluated first. For 
example, the logical expression 

( ( I .GT . (B+C) ) . AND. L) 

is evaluated in this order 



1. B+C 

2. I.GT.X 

3. Y.AND.L 



Call the result X 
Call the result Y 
Final operation 



The logical expression to which the logical operator .NOT. applies 
must be enclosed in parentheses if it contains two or more quantities. 
For example, assume that the values of the logical variables A and B are 
• FALSE, and .TRUE., respectively. Then these two expressions are not 
equivalent 

.NOT. (A. OR. B) 
. NOT. A. OR. B 

In the first expression, A. OR. B is evaluated first. The result is 
.TRUE.? but .NOT. (.TRUE.) implies .FALSE.. Therefore, the value of the 
first expression is .FALSE.. 

In the second expression, .NOT. A is evaluated first. The result is 
.TRUE.; but .TRUE. .OR. B implies .TRUE.. Therefore, the value of the 
second expression is .TRUE.. 
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ARITHMETIC AND LOGICAL ASSIGNMENT STATEMENT 



If 


J. 


W 




A # 


B, 


C r 


D 


h 








G„ 


H 







r ._. _ 1 

I General Form | 
|_„ . . 

a = b 

where a is any subscripted or nonsubscri pted variable 

b is any arithmetic or logical expression 

Note: a roust be a logical variable if , and only if, b is a logic- 

al expression. 

The FORTRAN arithmetic and logical assignment statement closely 
resembles a conventional algebraic equation; however, the equal sign of 
the FORTRAN arithmetic statement specifies replacement rather than equi- 
valence* That is f the expression to the right of the equal sign is eva- 
luated, and the resulting value replaces the current value of the vari- 
able to the left of the equal sign. 

Assume that the type of the following variables has been specified 

Varia ble Names Type Length Specification 

Integer variables **,<*, 2 

Real variables 4,4,8,8 

Complex variable 8 

Logical variables 4,4 

Then, the following examples illustrate valid arithmetic statements 
using constants, variables, and subscripted variables of different 

types « 

Statements Description 

A = B The value of A is replaced by the current value of B. 

*i = B The value of B is converted to an integer value and 

the least significant part replaces the value of W. 

A = I The value of I is converted to a real value and this 

result replaces the value of A. 

1=1+1 The value of I is replaced by the value of I + 1. 

E = I**J+D I is raised to the power J and the result is con- 

verted to a real value, to which the value of D is 
added. This result replaces the real part of the 
complex variable E. The imaginary part of the com- 
plex variable is set to zero. 

A = C*D The most significant part of the product of C and D 

replaces the value of A. 

G = .TRUE. The value of G is replaced by the logical constant 

.TRUE.. 

II = .NOT.G If G is .TRUE., the value of H is replaced by the 

logical constant .FALSE.. If G is .FALSE., the value 
of H is replaced by the logical constant .TRUE.. 
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G - 3..GT.I The value of I is converted to a real value? if the 

real constant 3. is greater than this result, the 
logical constant .TRUE, replaces the value of G. If 
3. is not greater than I, the logical constant 
• FALSE. replaces the value of G. 

E = CI. 0,2.0) The value of the complex variable E is replaced by 

the complex constant (1.0,2.0). Note that the state- 
ment E = (A,B), where A and B are real variables, is 
invalid. 

A = E The real part of the complex variable E replaces the 

value of A. 

E = a The value of A replaces the value of the real part of 

the complex variable Ej the imaginary part is set 
equal to zero. 
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CONTROL STATEMENTS 



Normally, FORTRAN statements are executed sequentially? that is, 
after one statement has been executed, the statement immediately follow- 
ing it will be executed. This section discusses the statements that may 
be used to alter and control the normal sequence of statement execution 
in the program. 

GO TO STATEMENTS 

The GO TO statements cause control to be transferred to the statement 
specified Dy a statement number. The three GO TO statements are: 
unconditional GO TO, computed GO TO, and assigied GO TO. Every time the 
same unconditional GO TO statement is executed, a transfer to the same 
specified statement is made. However, the computed and assigned GO TO 
statements cause control to be transferred to one of several statements, 
depending upon the current value of a particular variable. 

Unconditional GO TO Statement 



j General Form | 

| GO TO X | 

I I 

| where X is an executable statement number I 

This GO TO statement causes control to be transferred to the state- 
ment specified by the statement number. Every subsequent execution of 
this GO TO statement results in a transfer to that same statement. 

Example : 

50 GO TO 25 
10 A = B * C 



25 C = E**2 



Explanation ; Every time statement 50 is executed, control is trans- 
ferred to statement 25. 

Computed GO TO Statement 

| General Form | 

GO TO (x ± , x 2 , x 3 , ...,x n ), i 

where x ir x a ...,xn, are executable statement numbers 

i is a nonsubscripted integer variable in the range: 
1 > i > n 
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This statement causes control to be transferred to the statement num- 



bered x ±9 x 



at *3» • 



, or xn, depending on whether the current value of i 



is 1, 2, 3,..., or n f respectively. If the value of i is outside the 
allowable range, the next statement is executed. 



Example: 



GO TO (25, 10, 50, 7), ITEM 



50 A = 



C = 



25 L = 



B+C 



E**2*A 



C.GT.D.AND.F.LE.G 



10 B = 21.3E02 

Explanation : If the value of the integer variable ITEM is 1, statement 
25 will be executed next. If ITEM is equal to 2, statement 10 will be 
executed next, and so on. 



ASSIGN and Assigned GO TO Statements 



General Form 



1 

I 
4 



ASSIGN i TO m 

GO TO ra, (Xi,x a ,x 3l . . . ,xn) 

where i is an executable statement number 

*k*2£a»x 3 ,. • • *3£n are executable statement numbers 

m is a nonsubscripted integer variable of length <*, to which 
is assigned one of these statement numbers: Xi,x 2 ,x 3 ,...x n , 









The assigned GO TO statement causes control to be transferred to 
the statement numbered x JL ,x a ,x 3 ,. . . ,or x m depending on whether the 
current assignment of m is x if x a ,x 3 ,. . . ,or x n , For example, in the 
statement 

GO TO N, (10, 25, 8) 

if the current assignment of the integer variable N is statement 8, 
that statement is executed next. If the current assignment of N is 
statement 10, that statement is executed next. Similarly, if N is 
assigned statement number 25, that statement is executed next. 

The current assignment of the integer variable m is determined by 
the last ASSIGN statement executed. Oily an ASSIGN statement may be 
used to initialize or change the value of the integer variable m. The 
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v<»] uv of the Integer variable m is not the integer statement number; 
ASSIGN 10 TO I is not the same as 1=10. 



ASSIGN 50 TO NUMBER 

10 GO TO NUMBER, C35 f 50 # 25, 12, .18) 



\t) A = B + C 



,; M , .^.l J r J,t -_±.Vf l: Statement 50 is executed immediately after statement 10, 



ASSIGN 10 TO ITEM 



13 GO TO ITEM, (8, 12 r 2 5 , 50, 10) 



A « B + c 



b - C + D 

ASSIGN 25 TO ITEM 
GO TO 13 



C ~ l'**2 



? x L ) ]. , J. ri ^-A9i} : Th e first time statement 13 is executed r control is 

*- r.iiv.iter ced to statement 10. On the second execution of statement 13 # 

control it? transferred to statement 25. 

A r i *_ J 1 1 i'e lie IF St a t em ent 



j General Form | 

j IF (a) x rL# x 2# x 3 

1 

> where a is an arithmetic expression which is not complex 

* ^i^a^-i are statement numbers 

This statement causes control to be transferred to the statement num- 
s«.-r:ed x JL #x a ,x 3 when the value of the arithmetic expression a is less 
than, equal to, or greater than zero, respectively* The first execut- 
,-iblH statement following the arithmetic IF statement should have a sta- 
G«if",ui n.imber; otherwise, it can never be referred to or executed. 
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Example s 



IF (A(J,K)**3-B)10, <4, 30 



4 D = B ♦ C 



30 C = 



10 E = 



D**2 



(F*B)/D*1 



Explanation ; If the value of the expression (A(J,K>**3-B) is negative, 
statement 10 is executed next. If the value of the expression is zero, 
statement *4 is executed next. If the value of the expression is posi- 
tive, statement 30 is executed next. 

Logical IF Statement 






| General Form 

, _ _ _ „„ „„ _ 

IF(a)s 

where a is any logical expression 

s is any statement except a specification statement, DO sta- 
tement, or another logical IF statement 



The logical IF statement is used to evaluate the logical expression a 
and to execute or skip statement s, depending on whether the value of 
the expression is .TRUE, or .FALSE., respectively. 

Example 1 ; 



5 IF(A.LE.O.O) GO TO 25 

10 C = D ♦ E 

15 IF(A.EQ.B) ANSWER = 2.0*A/C 

20 F = G/H 



25 W - X**Z 



Explanation : In statement 5 , if the value of the expression is .TRUE. 
(i.e., A is less than or equal to 0.0), the statement GO TO 25 is 
executed next, and control is passed to statement 25. If the value of 
the expression is .FALSE. (i.e., A is greater than 0.0), the statement 
GO TO 25 is ignored, and control is passed to statement 10. 

In statement 15, if the value of the expression is .TRUE. (i.e., A 
is equal to B) , the value of ANSWER is replaced by the value of the 
expression (2.0*A/C), and statement 20 is executed. If the value of the 
expression is .FALSE. (i.e., A is not equal to B>, the value of ANSWER 
remains unchanged, and statement 20 is executed next. 
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Example 2 : Assume that P and Q are logical variables. 



5 IF(P.OR..NOT.Q)A=B 
10 C = B**2 



Explanation : In statement 5, if the value of the expression is .TRUE., 
the statement A=B is executed next and control continues to statement 
10. If the value of the expression is .FALSE., the statement A=B is 
skipped and statement 10 is executed. 

DO Statement 



r ^ — .„__. . . . . ^ 

1 General Form | 

| End of DO Initial Test 

| range variable value value Increment 



DO 



!i!a* 



( where x is the statement number of an executable statement that 
j follows the DO statement 

I 

| i is a nonsubscripted integer variable 

I 

1 5)i.» Ea' 5!3» are either unsigned integer constants greater 

| than zero or unsigned nonsubscripted integer variables whose 

1 values are greater than zero. The sum m 2 +m 3 +l must not 

| exceed the size of virtual storage. *m 3 , is optional; if it 

| is omitted, its value is assumed to be 1. In this case, the 

| preceding comma must also be omitted. > 

The DO statement is a command to execute repeatedly the statements 
that follow, up to and including the statement numbered x. The range of 
a DO is that set of statements that will be executed repeatedly; i.e., 
it is the sequence of consecutive statements immediately following the 
DO statemment. The first time the statements in the range of the DO are 
executed, i is initialized to the value m+; each succeeding time i is 
increased by the value m>. When, at the end of an iteration, i is equal 
to the highest value that does not exceed m< , control passes to the sta- 
tement following the statement numbered x. Thus, the number of times 
the statements in the range of the DO is executed is given by the 
expression : ° 

r i 

I — — — | +1 

1 5?3 I 
L J 

where the brackets represent the largest integral value not exceeding 
the value of the expression. If m a is less than m ± , the statements in 
the range of the DO are executed once. Upon completion of the DO, the 

DO variable is undefined. 

There are several ways in which looping (repetitively executing the 

same statements) may be accomplished when using the FORTRAN language. 
For example, assume that a manufacturer carries 1000 different machine 
parts in stock. Periodically, he may find it necessary to compute the 
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amount of each different part that is presently available. This amount 
may be calculated by subtracting the number of each item used f OUT(I), 
from the previous stock on hand, STOCKCIK 



Example 


1: 


5 


1=0 " 


10 


1=1*1 


25 


STOCKCD^STOCKCI)- OUTCI) 


15 


IFCI-1000) 10 f 30 f 30 


30 


A=B*C 



Explanation ; The three statements (5 f 10, and 15) required to control 
the loop could be replaced by a single DO statement , as shown in Example 
2. 

Example 2 i 



DO 25 I = l f 1000 
25 STOCK CI) = STOCKCI)-OOT(I) 
30 A=B*C 



Explanation : The DO variable, I f is set to the initial value of 1. 
Before the second execution of statement 25, I is increased by the 
increment 1 and statement 25 is again executed. After 1000 executions 
of the DO loop, I equals 1000. Since I is now equal to the highest 
value that does not exceed the test value, 1000, control passes out of 
the DO loop, and statement 30 is executed next. Note that the DO vari- 
able I is now undefined; its value is not necessarily 1000 or 1001. 

Example 3 : 



DO 25 1=1, 10, 2 
15 J=I*K 

25 ARRAY (J) = BRAY (J) 
30 A=B*C 



Explanation : Statement 25 is the end of the range of the DO loop. The 
DO variable, I, is set to the initial value of 1. Before the second 
execution of the DO loop, I is increased by the increment 2, and state- 
ments 15 and 25 are executed a second time. After the fifth execution 
of the DO loop, I equals 9. Since I is now equal to the highest value 
that does not exceed the test value, 10, control passes out of the DO 
loop, and statement 30 is executed next. Note that the DO variable I is 
now undefined; its value is not necessarily 9 or 11. 
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Programming Considerations in Using a DO Loop 



1. The indexing parameters of a DO statment (i, m x , m a , m 3 ) may not be 
changed by a statement within the range of the DO loop, or by any 
subprograms that are called within the range of a DO loop. 

2. A DO statement may contain other DO statements within its range. 
Ail statements in the range of an inner DO must be in the range of 
the outer DO. A set of DO statements satisfying this rule is 
called a nest of DOs. 

Example 1 ; 

DO 50 I = 1, U 

ACI) = BCI)**2 

DO 50 J=l, 5 

50 CCJ+1) = A(I> 

Example 2 : 

DO 10 INDEX = L, M 

N = INDEX ♦ K 

DO 15 J = 1, 100, 2 

15 TABLE (J) = SUM(J,N)-1 

10 BCN) = ACN) 



Range of 
inner DO 



Range of 
inner DO 



Range of 
outer DO 



Range of 
outer DO 



3. A transfer out of the range of any DO loop is permissible at any 
time - 

** . If, and only if, a transfer is made from the range of an innermost 
DO loop, transfer back into that loop is allowed, provided none of 
the indexing parameters (i,m 1 ,m 2 ,m 3 ) are changed outside the range. 
A transfer back into the range of any other DO within a nest of DOs 
is not permitted. 



Example : 



DO 



DO 



DO 







-► 2 



-> 



— V' 



-> 



Explanation ; The transfers specified in the example by the numbers 
1, 2, and 3 are permissible; those specified by **, 5, and 6 are 
not. 
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5, The indexing parameters (i f m^ ,m af ni 3 ) may be changed by statements 
outside the range of the DO statement only if no transfer is made 
back into the range of the DO statement using those parameters* 

6. The last statement in the range of a DO loop (statement x) may not 
be a GO TO f arithmetic IF, PAUSE, STOP, RETURN, or another DO sta- 
tement. Also, the last statement may not be a logical IF statement 
containing any of those statements, 

CONTINUE Statement 

I General Form 1 

| CONTINUE | 

CONTINUE is a dummy statement which may be placed anywhere in the 
source program without affecting the sequence of execution, it may be 
used as the last statement in the range of a DO statement to avoid end- 
ing the DO loop with any of the statements that are not permitted as the 
last statement in the range of a DO. 

Example 1 : 



DO 30 1=1, 20 
7 IF (ACl)-Btl)) 5,30,30 
5 ACI) =A(I) +1-0 

B(I) - BCD -2.0 

GO TO 7 
30 CONTINUE 
40 C = A(3) * BC7> 



Explanation : The CONTINUE statement is used as the last statement in 
the range of the DO statement, to avoid ending the DO loop with the sta- 
tement GO TO 7, 

Example 2 ; 



DO 30 1=1,20 

IF(A(I>-B(I>>5, 40,40 
5 A(I) = CCI) 

GO TO 30 
40 ACI) = 0,0 
30 CONTINUE 



Explanation : The CONTINUE statement provides a branch point that 
enables the programmer to bypass the execution of statement 40, 
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PAUSE Statement 

i 1 

I General Form | 

PAUSE 

PAUSE n 

PAUSE ' message ' 

where n is an unsigned 1-through 5- digit integer constant 

message is any literal constant 



A PAUSE statement executed in a program results in a message being 
written as follows: 

PAUSE Statement Resulting Message 

PAUSE PAUSE 00000 

PAUSE n PAUSE 1-5 digit integer 

PAUSE 'message' PAUSE text of message 

In nonconversational mode, the message is written on the standard 
system output data set and the pause is ignored, so the program con- 
tinues execution at the next executable statement. In converstaional 
mode of execution, the pause message is written at the user's terminal 
and the program waits until the user resumes execution via the TSS com- 
mand system, 

STOP Statement 

| General Form j 

V J 

j STOP | 

| STOP n 1 

I I 

| where n is an unsigned 1- through 5-digit integer constant 1 

L . . i 



This statement terminates the execution of the object program; mes- 
sage will be displayed as follows: 

STOP Statemen t Message 

STOP STOP 

STOP n STOP 1-5 digit integer 

END Statement 

r 1 

| General Form I 

I END 



The END statement is a nonexecutable statement that defines the end 
of a source program or source subprogram for the compiler. Physically, 
it must be the last statement of each program or subprogram. 

The END statement must be contained on a single line. 



-j 
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INPUT/OUTPUT STATEMENTS 



The input/output statements enable a user to transfer data, belonging 
to a named collection of data, between input/output devices and internal 
storage. The named collection of data, called a data set, is a con- 
tinuous string of data that may be divided into FORTRAN records. 

A data set is referred to by an integer constant or integer variable, 
called the data set reference number. 

There are five I/O statements: READ, WRITE, END FILE, REWIND, and 
BACKSPACE. The READ and WRITE statements cause transfer of records from 
and to data sets and internal storage. The END FILE statement defines 
the end of a data set; the REWIND and BACKSPACE statements control the 
positioning of data sets. 

In addition to these five statements, the FORMAT and NAMELIST state- 
ments, although they are not I/O statements, are used with certain forms 
of the READ and WRITE statements. The FORMAT statement specifies the 
form in which the data is to be transmitted! the NAMELIST statement spe- 
cifies a list of variables or array names to be used in an input/ output 
operation. Also, both statements allow the user to divide a data set 
into FORTRAN records. 

Even though the I/O statements are device independent, the source or 
the destination of the data being transferred influences the specifica- 
tion of the records and data formats. Therefore, subsequent examples 
are in terms of card input and print-line output, unless otherwise 

noted. 



READ STATEMENT 



| General Form 



READ (a, b, END=C, ERR=d) list 

where a is an unsigned integer constant or an integer variable of 
length ** that represents a data set reference number. 

b is either the statement number or array name of the FORMAT 
statement describing the data being read, or a NAMELIST 

name. 

c is the statement number to which transfer is made upon 
encountering the end of the data set. 

d is the statement number of the statement to which transfer 
is made upon encountering an error condition in data 

transfer . 

list is a series of variable or array names, separated by 
commas , which may be indexed and incremented? they specify 
the number of items to be read and the storage locations 
into which the data is placed. 



The READ statement may take many different forms. For example, the 
parameters END=c and ERR=d are optional and, therefore, may or ma y not 
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■ un^Mi jr. '. HPAA .- <■ ih'i'rut . r* ; rt h* 1 r more tf the parameter l_i£t # the para- 
">*<?, r _ . r - .,. _ ,, , , , . i , , ,, f , ; J < , .,^ 

'" ; 0'O • '-• : - •* •'• Oo , >•> i -or*- o oo, FNb=c or ERR=d is used after the a 
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.cad o< t ,!;0 soioi- . f>o»'u',e.rt'i.fi:; ! he end of the data set. if an END 
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term ;. it ■ -♦ . »-. ^.-u \r,vi -.f the object program, 
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no to, : , ■ .■ '<':'■>. \ ?•• ■« i •, o- i -"e : !io associated with the record in error. 
Ac i \ • o , . ',...• . , :=,, oo> p • • op am as to which input record or reco- 
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lata o\ 4 ~ ? I i A-- !":;■. AD li^t. "if an ERR parameter is not specified in a 
;OAn t , : j »•-. ,f^o , ;U .. v ? tJl t <>rm mater, execution of the object program. 

"i'\ *' h ) : --* .- .i5-- , ,"' ioTfit; -'A; the REftD statement are 



o.'A5Aa}!, o 

T v ,* par aoe i er y>. f^'! v ~,'; and FRR-ri nay be used, in the combination 

;cr • A< d a At o op A" *• : ■■ n of these three forms. 



Ti> ' ; v';: ; ,• • •;•>->.; "■ • '■ ioid oat a front the data set associated with a 
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>;Tr.-/ natr.cs t "^reives a 'NAMELIST name by use of a 
'l t r programmer need only use the NAMELIST name in 
•:f'-.'- i- " r^feu'iw^ that list thereafter in the 
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i :A*. ..':': I : . !"/'- "'-», f :>, . ,c /y/A ,e . . . , r f/2/g # h r . . • ,i 

' " I 

1 ..Ai»'. .-- - , ( ,o> j v., . , .i/e NAMELIST names | 

I 

| .1, t, '",o^ « • - are variable or array names | 

Tn-? tA^lA;^i?s i.\s 1 ^-s riorAiy to defining and using a NAMELIST name; 

1, A N/ r ,u o.j A*, ( r, ^^ oo>.sist.5i of from 1 through 6 alphameric characters, 
■, .,•.. •• ,■. i o ' ;• a'A; A : , alphabetic. 

2. -\ f'i/A-A: ,i,j : '' ^auu, i^» e^clo:>ed in slashes. The list of variable or 
oro-o n ; '.."! '.•? ' ' lioiq.o.oj to a NAMELIST name ends with a new NAMELIST 

.•o- •■!!>• o^-v • o> c - ».--"o-h' ; o; or with the end of the NAMELIST statement. 



3. A variable name or an array name may belong to one or more NAMELIST 
names . 

U • A NAMELIST name may be defined only once by its appearance in a 
NAMELIST statement and must be so defined before its use. the 
NAMELIST name may appear only in READ or WRITE statements in tiie 
program. 

5. A NAMELIST statement may appear anywhere in a FORTRAN program prior 
to its use in a READ/WRITE statement. 

6. Variable or array names appearing anywhere in a NAMELIST statement 
or NAMELIST name may not appear in a FUNCTION, SUBROUTINE, or ENTRY 
statement. 

Example : Assume that A, I, and L are array names. 



NAMELIST /NAM1/A,B,I,J,L/NAM2/A,C, J f K 



READ (5,NAM1) 



Explanation : The READ statement causes the record that contains the 
input data for the variables and arrays that belong to the NAMELIST 
name, NAMl, to be read from the data set associated with the data set 
reference number 5. 

Input Data 

When a READ statement refers to a NAMELIST name, input data in the 
form described below is read from the designated input data set. 

The first character in the record must always be blank. The second 
character of the first record of a group of data records to be read must 
be & (ampersand), immediately followed by the NAMELIST name. The NAME- 
LIST name must be followed by a blank and must not contain embedded 
blanks. This name is followed by any combination of data items 1 and 2 
below, separated by commas. (A comma after the last item is optional.) 
The end of a data group is signaled by £END. 

The form the data items may take is 

1 . Variable name = constant 

Hie variable name may be a subscripted variable name or a single 
variable name, subscripts must be integer constants. 

2. Array name = set of constants (separated by commas) 
The set of constants may be in the form 

k* constant 

where k is an unsigned integer called the repeat constant. It 
represents the number of successive elements in the array to be 
initialized by the specified constant. The number of constants 
must not be greater than the number of elements in the array. 
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Input constants may also be Hollerith (H format) or hexadecimal (Z 
format) data. The H format is used as in FORMAT statements. The repeat 
constant may not be used with the H- format option. The size of the 
character string should not exceed the size o f an element. To use the Z 
format, prefix the hexadecimal characters to be read with a "Z". 



Constants used in the data items may be integer, real, literal, com- 
plex, or logical data. If the constants are logical data, they may be 
in the form T or .TRUE. and F or .FALSE.. 

Any selected set of variable or array names belonging to the NAMELIST 
name appearing on the first record may be used as specified by items 1 
and 2 in the preceding text. Names that are made equivalent to these 
names may not be used unless they also belong to the NAMELIST name. 

The end of a group of data is signaled by the character string SEND, 
with no embedded blanks and all appearing in the same record. 

Blanks must not be embedded in a constant or repeat constant, but may 
be used freely elsewhere in a data record. The last item on each record 
that contains data items must be a constant followed by a comma. (The 
comma is optional on the record that precedes the record containing 
&END.) Trailing blanks after integers and exponents are treated as 
zeros . 

Example ; Assume that L is an array consisting of one subscript paramet- 
er ranging from 1 to 10. 

Column 2 
■ 
First data card: SNAMl I<2,3>=5, J=4, 



Last data card: A(3)=4.0, L=2, 3, 8* 4, SEND 

Explanation : If this data is input to be used with the NAMELIST and 
READ statements previously illustrated, the following actions take 
place: The first data card is read and examined to verify that the name 
is consistent with the NAMELIST name in the READ statement. If the name 
does not match the NAMELIST name, the next NAMELIST group is read. Vfoen 
the data card is read, the integer constants 5 and 4 are placed in 
1(2,3) and J, respectively; the real constant 4.0 is placed in A (3). 
Since L is an array not followed by a subscript, the entire array is 
filled with the succeeding constants. Therefore, the integer constants 
2 and 3 is placed in L(l) and L(2), respectively, and the integer con- 
stant 4 is placed in L(3), L( 4) ,. . . ,L( 10) . 

READ (a,b) List 

This form is used to read data from the data set associated with a 
into the storage locations specified by the variable names in the list . 
Tne list , used in conjunction with the specified FORMAT statement b (see 
"FORMAT statement"), determines the number of items (data) to be read, 
the locations, and the form the data will take in storage. 

Example 1 : Assume that the variables A, B, and C have been declared as 
integer variables. 



75 FORMAT (G10, G8 , G9) 
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READ (J, 75) A, B f C 



Explanation : The READ statement above causes input data from the data 
set associated with data set reference number J to be read into loca- 
tions A, B f and C, according to the FORMAT statement referenced (state- 
ment 75). That is f the first 10 characters of the record are read, con- 
verted to internal form, and stored into A, the next 8 characters into 
3, and the next 9 characters into C. 

Tne list can be omitted from the READ (a, b) list statement. In this 
case, a record is skipped or data is read from the data set associated 
with a, into the storage locations occupied by FORMAT statement b. 

Example 2 : 



98 FORMAT ( , HEADING•) 



READ (5,98) 



Explanation ; The statements above would cause the characters H, E, A, 
D, I, N, and G, in storage, to be replaced by the next 7 characters in 
the data set associated with data set reference number 5. 

Example 3 : 



98 FORMAT (G10 , < HEADING 1 ) 



READ (5,98) 



Explanation ; The statements above would cause the next record in the 

data set associated with data set reference number 5 to be skipped. No 

data is transferred into internal storage because there is no list item 
that corresponds with format code G10. 

READ (a) List 

The form READ (a) list of the READ statement causes binary data 
(internal form) to be read from the data set associated with a into the 
storage locations specified by the variable names in the list . Since 
the input data is always in internal form, a FORMAT statement is not 
required. This statement is used to retrieve the data written by a 
WRITE (a) list statement. 

Example; 

READ (5) A, B, C 
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Explanation : This statement causes the binary data from the data set 
associated with data set reference number 5 to be read into the storage 
locations specified by the variable names A, P, and C. 

r ^ ie li st may be omitted from the READ (a) list statement; in this 
case, a record is skipped. 

Example ; 

READ (5) 

Explanation : The statements above would cause the next record in the 
data set associated with data set reference number 5 to be skipped. No 
data is transferred into internal storage. 

Index ing I/O Lists 

Variables within an I/O list may be indexed and incremented in the 
same manner as those within a DO statement. These variables and their 
indexes must be included in parentheses. For example, suppose it is 
desired to read data into the first five positions of array A. This may 
be accomplished by using an indexed list: 

15 FORMAT CGI 0.3) 



READ (2 # 15> (A(I),I=1,5) 
This is equivalent to 
15 FORMAT (G10.3) 



DO 12 I = 1,5 
12 READ (2, 15) A( I> 

As with DO statements, a third indexing parameter may be used to spe- 
cify the amount by which the index is to be incremented at each itera- 
tion. Thus, 

READ (2,15) CACI), 1=1,10,2) 
causes transmission of values for AC1), A(3), A(5), A(7), and A(9). 
Furthermore, this notation may be nested. For example, the statement 

READ (2,15) (CC(I,J),D(I,J) , J=l ,3 ) , 1=1 ,4 ) 

would transmit data in this order: 

C(l,l), D(l,l), C(l,2), D(l,2), C(l,3), D(l,3) 

C(2,l), D(2,l), C(2,2), D(2,2), C(2,3), D(2,3) 

C(3,l), D(3,l), C(3,2), D(3,2), C(3,3), D(3,3) 

CU,1), D(4,l), C(4,2), D(4,2), C(4,3), D(4,3) 

Since J is the innermost index, it varies more rapidly than I. 

As another example, consider 

READ (2,25) I , (C ( J) , J=l ,1 ) 

The variable I is read first and its value then serves as an index to 
specify the number of data items to be read into array C. 
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If it is desired to read data into an entice arc, M it. l .v rrjt n»-eef-- 
sary to index that array in the I/O list. Fox. ^x.^-ie, us-.nm.e \ t^* the 
array A consists of one subscript parameter, varying in t h*- r:vi'-" ot 1 
to 10. Then this READ statement, referring ;.u F\>KMA ? ..UUi^f*/ asmbered 

READ (2,5) A 

would cause data to be rc?ad into h (1) f A ( 2) # .. . . ,A il i „ 

The indexing of I/O lists applies to WHITK ii.-, t,: r ,* - v-1 ,t- klad 

lists . 

Reading Form at Statements 

FORTRAN provides the facility for var iabj e W-MAI :o:.-i t ^ooi s oy 
allowing a FORMAT statement to be zead into an ,-irnn, hi i.'is-r.v.H' .?. na 

using the data in the array as the FORMAT s (*><-i i - • -:i # m.sls t^:r w'U'iihju^tu 
I/O statements « 

For example, the statements below result if., ft, E, and diniy C being 
read, converted, and stored according to the* format ,* peed f icat i ons read 

into the array FMT at object time. 

DIMENSION FMT (18) 
1 FORMAT (18A4) 

READ C5 # l> FMT 
READ C5 # FMT) A,B, (C( I ) , 1=1 9 5 ) 

1. The name of the variable forma*- speci t i <:<**-. ion ^ ,: - : a > ;.?•£ »■ ,; n .; 

DIMENSION statement, even if the array ;<i/,e a-, juj.' i. 



K 



The form of the format codes read into c he KC -o. ?,>i A *i 

must take the same form as a source program Fv.-MiVV st«* foment * 

except that the word FORMAT is omitted (see **'•;{ -e <*T>*;v.//r 

Statement 88 ) . 

If a format code read in at objee*: i: iroe a^M.ai is d«ni:i.N*< • q^.-» r* ophe s 
within a literal field that is defcined by apostrophes, i ? .in o a id be 
used for output only. If an object tiioe foima-. rode is to i->^ used 
for input and if it must contain a Litetrai. Lie id wirn an itii un-ii 
apostrophe, the H format code must be used tcr t.he literal t x«r id 
definition. 



WRITE STATEMENT 



*i 



| General Form j 

WRITE <a r b) list | 

I 

where a is an unsigned integer const a m~ or -\u -"ilt-ge; vitutlr- of [ 

length 4 that represents a data ser refe-^nce in*:i;»v?f, j 

i 
b is either the statement number cr -l^rn, name r* -'.*• FORMAT | 
statement describing the data ty^im:; w* ^ c -"n, oi ,:> NAME! X-ST j 

name • j 

I 

list, is a series of variable or .:r.rviay nvies, sop.ar .red cy f 

commas, which may be indexed and in oi ew/ ■ t ed ; i-nov specify | 

the number of items to be written and the ntor i.jc l 'K>,t ions j 
| from which 'the data is taken. j 
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The WRITE statement may take many different forms. For example, the 
list or the parameter b may be omitted. 

The three basic forms of the WRITE statement are 

WRITE(a,x) 

WRITE(a,b)list 

WRITE(a) list 

WRITE (a g x) 

This form is used to write data from the storage locations specified 
by the NAMELIST name x into the data set associated with a (see 
"READ (a, x)") . 

Example : 

WRITE(6,NAM1) 

Explanation ; 

Ihis statement causes all variable and array names (as well as their 
values) that belong to the NAMELIST name, NAM1, to be written on the 
data set associated with data set reference number 6. 

When a WRITE statement references a NAMELIST name 

1. All variables and arrays and their values belonging to the NAMELIST 
name will be written out, each according to its type. The complete 
array is written out by columns, 

2. The output data will be written such that 

a. The fields for the data will be large enough to contain all the 
significant digits; 

b. The output can be read by an input statement referencing the 
NAMELIST name. 

Example : Assume that A is a 3-by-3 array. 



NAMELIST/NAM1/A,B,I,D 
WRITE (8 r NAMl) 



Assuming that the output is punched on cards, the format would be: 

Output Card Column 2 

■ 

First SNAM1 

Second A=3.4, 4.5, 6.2, 25.1, 

Third 9.0, -15.2,-7.6, 0. 576Ebl2, 

Fourth 2. 717, B=3. 14,I=10,D=0. 378E-15, 

Fifth SEND 

WRITE (a,bJ List 

This form is used to write data in the data set associated with a 
from the locations in storage specified by the variable names in the 
list . The list , used in conjunction with the specified FORMAT statement 
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b, determines the number of items (data) to be written, the locations, 
and the form the data will take in the data set. 

Example 1 : 

75 FORMAT (G10, G8, G9) 



WRITE (J f 75) A f B # C 

Explanation ; The WRITE statement above causes output data to be written 
in the data set associated with the data set reference number J f from 
locations A f B f C, according to the FORMAT statement referred to (state- 
ment 75). (Format statements are described in a later section.) 

The list may be omitted from the WRITE (a f b) list statement. In this 
case, a blank record is inserted, or data is written in the data set 
associated with a from the locations in storage occupied by FORMAT sta- 
tement b. 

Example 2 : 

98 FORMAT ( - HEADING") 



WRITE (5,98) 

The statements above would cause a blank and the characters H, E, A, 
D, I, N, and G in storage to be written in the data set associated with 
data set reference number 5. 

Example 3 : 

98 FORMAT (G10, •HEADING") 



WRITE (5,98) 

Explanation ; The statements above would cause a blank record to be 
placed in the data set associated with data set reference number 5. No 
data is transferred into the data set because there is no list item that 
corresponds with format code G10. 

WRITE (a) List 

The WRITE (a) list form of the WRITE statement causes binary data 
(internal form) from the storage locations specified by the variable 
names in the list to be written in the data set associated with a. 
Since the output data is always in internal form, a FORMAT statement is 
not required. The READ (a) list statement is used to retrieve the data 
written by a WRITE (a) list statement. 

Example: 

WRITE (5)A f B, C 

Explanation : The statement causes the binary data from the locations 
specified by variable names A, B r and C to be written in the data set 
associated with data set reference number 5. 
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FORMAT STATEMENT 



1 General Form | 

| x FORMAT Cc iJP c a# • - . ,c n /c ± f ,c a f ,. . . #£n §/ - * •> 

i 

1 where x is a statement number CI through 5 digits) 

1 £a.*£a» • • • '£n an< ^ £i* *£%* *• • - *£n 8 are f° rmat codes which may 

| be delimited by one of the separators; comma, slash, or 

1 parenthesis; these codes specify the length, decimal point 

| (if any), and position of the data in the data set. 

I 

| The character /is used to separate FORTRAN records. 

i_„ .. __„_„„ .__. . ._ . . -_. , _. . j 

The FORMAT statement is used in conjunction with the READ and WRITE 
statements to specify the desired form of the data to be transmitted? 
the form is varied by the use of different format codes. 

The format codes are 

G -- to transfer integer, real, complex, or logical data 

I -- to transfer integer data 

F -- to transfer real data that does not contain a decimal exponent 

D -- to transfer real data that contains a D decimal exponent 

E — to transfer real data that contains an E decimal exponent 

L -- to transfer logical data 

Z -- to transfer hexadecimal data 

A -- to transfer alphameric data 

Literal -- to transfer a string of alphameric and special characters 

H — to transfer literal data 

X -- to either skip data when reading or insert blanks when writing 

T — to specify the position in a FORTRAN record where transfer of data 

is to start 

P --to specify a scale factor 

Kny number used in a FORMAT statement, except the statement number or a 
I iteral, must be less than or equal to 255. 

USE OF THE FORMAT STATEMENT : This section contains general information 

on the FORMAT statement. The points discussed below are illustrated by 

the examples that follow. 

1. FORMAT statements are nonexecutable and may be placed anywhere in 

the source program. 

2- A FORMAT statement may be used to define a FORTRAN record, as 

follows ; 

a. If no slashes or additional parentheses appear within a FORMAT 
statement, a FORTRAN record is defined by the beginning of the 
FORMAT statement Cleft parenthesis) to the end of the FORMAT 
statement (right parenthesis). Thus, a new record is read when 



the format control is initiated (left parenthesis); a new reco- 
rd is written when the format control is terminated (right 
parenthesis) . 



Example : 



x FORMAT ( , -— - , ) 



--corresponds to 1 
FORTRAN record 



If slashes appear within a FORMAT statement, FORTRAN records 
are defined as starting from the beginning of the FORMAT state- 
ment and proceeding to the first slash in the FORMAT statement, 
from one slash to the next succeeding slash, or from the last 
slash to the end of the FORMAT statement. Thus, a new record 
is read when the format control is initiated, and thereafter a 
record is read upon encountering a slash; a new record is writ- 
ten upon encountering a slash or when format control is 
terminated. 

Example : 

x FORMAT ( / — « / ) 

< > <-—> < — -> 

1 I 1 

I I I 

each corresponds to 1 FORTRAN record 

If more than one level of parentheses appears with a FORMAT 
statement, a record is defined by the beginning of the FORMAT 
statement and the end of the FORMAT statement. At this point, 
the definition of the FORTRAN record continues at the first- 
level left parenthesis that is closest to the right of the FOR- 
MAT statement and finishes at the end of the FORMAT statement. 



Example 1 ; 



X FORMAT ( — - 



1 
( — 



2 
C — 



21 
)) 





— ) 

— > 

— > 



each corresponds to 1 FORTRAN record 



Example 2 : 



Oil 110 
x FORMAT (- — ( — -) - — C — ) ) 



-> 
-> 



each corresponds to 1 FORTRAN record 
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When defining a FORTRAN record by a FORMAT statement, it is impor- 
tant to consider the original source (input) or ultimate destina- 
tion (output) of the record. For example, if a FORTRAN record is 
to be punched for output, the record should not be greater than 80 
characters. For input, the FORMAT statement should not define a 
FORTRAN record longer than the actual record in the data set. 

3. Blank output records may be introduced or input records may be 
skipped by using consecutive slashes (/) in a FORMAT statement. If 
there are n consecutive slashes at the beginning or end of a FORMAT 
statement, n input records are skipped or n blank records are 
inserted between output records, respectively. If n consecutive 
slashes appear anywhere else in a FORMAT statement, the number of 
records skipped or blank records inserted is n-1. 

4. Successive items in an I/O list are transmitted according to suc- 
cessive format codes in the FORMAT statement, until all items in 
the list are transmitted. If there are more items in the list than 
there are codes in the FORMAT statement, control transfers to the 
preceding left parenthesis of the FORMAT statement, and the same 
format codes are used again with the next record. If there are 
fewer items in the list, the remaining format codes are not used. 

5 . A format code may be repeated as many times as desired by preceding 
the format code with an unsigned integer constant. 

(» . h limited parenthetical expression is permitted to enable repeti- 
tion of data fields according to certain format codes within a 
longer FORMAT statement. Two levels of parentheses, in addition to 
the parentheses required by the FORMAT statement, are permitted. 
The second level of parentheses facilitates the transmission of 
complex quantities. 

7 . When transferring data on input or output, the type of format code 
used, type of data, and type of variables in the I/O list should 
correspond. 

8. In the examples below, the output is shown as a printed line. A 
carriage control character - x' (see "Carriage Control w ) is speci- 
fied in the FORMAT statement but does not appear in the first print 
position of the print line. This carriaqe control character 
appears as the first character of the output record on any I/O 
medium except the printed line. 

G Fo r mat Code 

r —. — „_. . _„ — . — . . 1 

| General Form | 

l_„ _„„„_„ .„_ . __ . ___ 4 

aGw.s 

where a is optional and is an unsigned integer constant, less than 
or equal to 255, used to denote the number of times the same 
format code is repetitively referenced 

w is an unsigned integer constant, less than or equal to 
255, specifying the total field length 

s is an unsigned integer constant specifying the number of 

significant digits 

L . . — . . - 

The G format code is a generalized code, in that it may be used to 
determine the desired form of data, whether integer, real, complex, or 
logical . 
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The .s portion may be omitted when transferring integer or logical 
data; if present, it is ignored. When real data is transferred, the w 
portion of the G format code includes four positions for a decimal 
exponent field. 

If the real data, say n, is in the range 0.1<n<10**s (where s is the 
s portion of the format code Gw. s) , this exponent field is blank. 
Otherwise, the real data is transferred with an E or D decimal exponent, 
depending on the length specification (either four or eight storage 
locations, respectively) of the real data* 

If insufficient positions are reserved by s, the number is rounded to 
s digits for output, and truncated to s digits for input. If excessive 
positions are reserved by s, zeros are filled in on the right. 

For simplification, the following examples deal with the printed 
line; however, the concepts developed apply to all input/output media. 

Example 1 : Assume that the variables A, B, C, and D are real, with 
values of 292.7041, 82.43441, 136.7632, and .8081945, respectively. 

1 FORMAT Cx' ,G12.4,G12.5,G12.4,G12.7) 

2 FORMAT Cx s ,^13.4,G13. 5,G13. 4) 

3 FORMAT ("x f ,G13.4> 



WRITE (5, n) A, B, C, D 



Explanation ; 



If n had been specified as 1, the printed output would be (b repre- 
sents a blank) 

Print position 1 Print position 48 

■ ■ 

bbb29 2.7bbbbbb82.434bbbbbbbl36.8bbbb.8 0819 45bbbb 

If n had been specified as 2, the printed output would be 

Print position 1 Print position 39 

bbbb292.7bbbbbbb82.434bbbbbbbbl36. 8bbbb Line 1 

bbbO.80 82bbbb Line 2 

It can be seen that by increasing the field width reserved (w), 
blanks are inserted. 

If n had been specified as 3, the printed output would be 

Print position 1 

bbbb292.7bbbb Line 1 

bbbb8 2.43bbbb Line 2 

bbbbl36.8bbbb Line 3 

bbb0.8082bbbb Line 4 

The same format code was used for each variable in the list. 
Each repetition of the same format code caused a new line to be 
printed. 
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Example 2: Assume that the variables I # J, K r and L are integer, with 
values of 292, 443428, 490B081 9 and 40018 # respective ly. 

1 FORMAT C * x 9 ,G10,2G7,G5) 

2 FORMAT Cx f ,G6) 

3 FORMAT C f x f , 4G10) 



WRITE (5 # li) I, J, K, L 

Explanat i o o : 

a. If n had been specified as 1, the printed output would be 
Print pos it ion 1 Fr in t po s i t i on 29 

• m 

bbbbbbb292b443428490808140018 Line 1 

The same results would be achieved, if FORMAT statement 1 had been 

wr it t eti as 

FORMAT ("x'jGIO, G7 # G7, G5) 

Note that the -s portion of the G format may be omitted when trans- 
mitting i n te g e r d a ta . 

o. If n bad been specified as 2, the printed output would be 

bbb?.92 Line 1 

443428 Line 2 

**♦♦♦+ Line 3 

b40018 Line 4 

NwU: f »i*it i he second rormat code G6 is an incorrect specification 
tor tne third variable K, i.e., 4908081. The field will be filled 
with asteri sks. 

c. If ii had been specified db 4, the printed output would be 

Pi ii 1 1. i K) si tion 1 P r i nt po si ti on 4 

bbbbbbb292bbbb4434 28bbb49080 81bbbbb40018 Line 1 

From the above example, ii onn be seeii that increasing the field 
width w improves readability. 



Kxarapi e i: Assume 



" J " 


Integer 


Lenqth 
2 


Value 
29 2 


A 


Real 


n 


4 71.93 


B 


Real 


a 


81.91 


D 


Heal 


8 


6.9310072 


C 


Complex 


8 


C2.1 # 3«7) 


L 


L<>7 i cal 


1 


.TRUE. 



1 FORMAT C ' x . * ,G 3 , 2G 9 . 2 ,G1 3.7, 2G8 . 2 , G3 ) 

2 FORMAT C f x* , G 3/ • x • , 2G10. 2/ ' x f # G9. 1 / f x s , 2G 8. 2, G3 ) 

3 FORMAT (//' x* ,G3, 2G9. 2//'x f ,G13. 7 , 2G 8 . 2 # G3///) 



a«4 



WRITE (5,n) I # A # B,D,C f L 



Explanation s 

a. If n has been specified as 1, the printed output would be 

Print position 1 Print position 53 

■ m 

292b0.47Eb03bb82.bbbbb6.931007bbbbb2.1bbbbb3.7bbbbbbT 

When complex data is being transmitted f two format codes are 
required* The real and imaginary parts are each treated as separ- 
ate real numbers, and the parentheses and comma are not printed as 
part of the output. 

b. If n has been specified as 2 f the printed output would be 

Print position 1 

■ 

29 2 Line 1 

bb(K«f7Eb03bbb82.bbbb Line 2 

bbb7«bbbb Line 3 

b2.1bbbbb3«7bbbbbbT Line <* 

The use of the slash (/) to separate two format codes causes the 
data not yet printed to be printed on a new line. If the output 
data is to be punched on cards* the slash specifies that the fol- 
lowing data will be punched on another card. 

c. If n has been specified as 3, the printed output would be 

Print position 1 

■ 

(blank line) 

(blank line) 

292b0.47Eb03bb82.bbbb 

(blank line) 

b6„931007bbbbb2. Ibbbbb3, TbbbbbbT 

(blank line) 

(blank line) 

(blank line) 

Note that the two consecutive slashes appearing at the beginning 
and the three at the end of the series of format codes cause blank 
lines to be inserted as shown* However, n consecutive slashes 
appearing elsewhere in a FORMAT statement cause the insertion of 
n-1 blank lines f as shown in line 4. 

The principles illustrated in the previous output examples also apply 
when using the READ statement on input- Also, there are further consi- 
derations when using the FORMAT statement on input or output. 

1. When reading real input data with a G format code, a decimal point 
must be included. 

2. The use of additional parentheses (up to two levels) within a FOR- 
MAT statement is permitted to enable the user to repeat the same 
format code when tx nsnrLtting data* For example, 
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Line 


1 


Line 


2 


Line 


3 


Line 


4 


Line 


5 


Line 


6 


Line 


7 


Line 


8 



10 FORMAT C2(G10.6,G7. 1) ,G*0 

is equivalent to 

10 FORMAT (G10.6, G7.1, G10.6, G7.1, G4 ) 

3. If a multiline listing is desired, with the first two lines to be 
printed according to a special format and all remaining lines 
according to another format, the last format code in the statement 
should be enclosed in a second pair of parentheses. For example, 

FORMAT Cx' ,G2,2G3.1/ , x , ,G10.8/( , x' ,3G5. . 1)) 

If more data items are to be transmitted after the format codes 
have been completely used, the format repeats from the last left 
parenthesis. Thus, the printed output would take the form 

G2,G3. 1,33.1 
G10.8 

G5.1,G5.1,G5.1 
G5.1,G5.1,G5.1 



As another example, consider the statement 

FORMAT ( 1 x , ,G2/2( i x < ,G3 ,G6. 1) ,G9. 7 ) 

If there are 13 data items to be transmitted, the printed output on 
a WRITE statement would take the form 

G2 

G3,G6.1, , x i ,G3,G6.1,G9.7 

G3,G6.1, f x' ,G3,G6.1,G9.7 

G3,G6.1 

Numeric Format Codes (I,F,E,D, and Z) 

Five types of format codes are available for the transfer of numeric 
data. These are specified in this form 



f 1 

[ General Form | 

alw 

aFw.d 

aEw.d 

aDw.d 

aZw 

where a is optional and is an unsigned integer constant, less than 
or equal to 255, used to denote the number of times the same 
format code is repetitively referenced 

I, F, E, D, and Z are format codes 

w is an unsigned integer constant less than or equal to 255, 
specifying the total field length of the data 

d is an unsigned integer constant specifying the number of 
decimal places to the right of the decimal point, i.e., the 
fractional portion 



46 



For purposes of simplification, the following description of format 
codes deals with the printed line. The concepts developed apply to all 
input/output media. 

I Format Code 

The I format code is used to transmit integer data. 

If the number of characters to be transmitted is greater than w, on 
input, the excess rightmost characters are lost? on output, the entire 
field, w characters, will be filled with asterisks. 

If the number of characters is less than w, on input, leading blanks 
are not significant? embedded and trailing blanks are treated as zeros. 
On output, the leftmost nositions are filled with blanks. 

If the quantity is negative, the position preceding the leftmost 
digit contains a minus sign. In this case, an additional position 
should be specified in w for the minus sign. If w is such that no space 
exists for the minus sign, the entire field, w characters, will be 
filled with asterisks. 

The following examples show how each of the quantities on the left is 
printed according to the format code 13 Cb represents a blank). 

Internal Value Printed Value 

721 721 

-7 21 *** (incorrect because of insufficient 

specification) 
-12 -12 

568114 *** (incorrect because of insufficient 

specification) 
bbO 

-5 b-5 

9 bb9 

F Format Code 

For F format codes which are used in conjunction with the transfer of 
real data, w is the total field length reserved, and d is the number of 
places to the right of the decimal point (the fractional portion) . This 
differs from the G format code, where the number of significant digits 
is specified. The total field length reserved must include sufficient 
positions for a minus sign (if any) and a decimal point. The sign, if 
negative, is printed. 

If insufficient positions are reserved by d, the fractional portion 
is rounded to the dth position. If excessive positions are reserved by 
d, zeros are filled in on the right. The integer portion of the number 
is handled in the same manner as numbers transmitted by the I format 

code. 

The following examples show how each of the quantities on the left is 
printed according to the format code F5. 2. 

Internal Value Printed Value 

12.17 12.17 

-41.16 ***** (incorrect, insufficient 

s pe ci f i cat ion ) 
-.2 -0.20 

7.354 2 b7.35 (last two digits of accuracy lost; 

insufficient specification) 
-1. -1.00 

9.03 b9.0? 

187 .64 ***** (incorrect; insufficient 

specification) 

Input/Output Statements 47 



D and E Format C&di ss 

Hie D and E format codes are used in conjunction with the transferral 
of real data that contains a D or E decimal exponent, respectively, A D 
format code indicates a field length of 8| an E code indicates a length 
of 4. For D and E format codes, the fractional portion is again indi- 
cated by d. The w includes field d, spaces for a sign, the decimal 
point, plus four spaces for the exponent* 

For output, space for at least one digit preceding the decimal point 
should be reserved., In general, w should be at least equal to d+7. If 
insufficient positions for d are supplied, the fraction is rounded to 
the dth position. If excessive positions are supplied, zeros are added* 

The exponent is the power of 10 by which the number must be multip- 
lied to obtain its true value. The exponent is written with a D or an 
E, followed by a space for the sign and two spaces for the exponent 
(maximum is 75). 

The following examples show how each of the quantities on the left is 
printed according to the format codes ( D10. 3/E10. 3) . 

Internal Value Printed Value 

238. b0„238Db03 

-.002 -G.200E-02 

.00 00 000000 4 bO .4000- 10 

-21.0057 -»0.210Eb02 (last three digits of accuracy 

lost; insufficient field width) 

When reading input data, the start of the exponent field must be 
marked by an E or, if that is omitted, by a ♦ or - sign (not a blank) . 
Thus, E2, E+2, + 2, +02, E02, and E + 02 all have the same effect and are 

permissible decimal exponents for input. 

Numbers for E, D, and F format codes need not have their decimal 
point punched. If it is not present, the decimal point is supplied by 
the d portion of the format code* If it is present in the card, its 
position overrides the position indicated by the d portion of the format 

code. 

Z Format Code 



| General Form | 

I ^ z ^ 

1 

| where a is optional and is an unsigned integer constant, less than 

I or equal to 255, used to denote the number of times the same 

1 format code is repetitively referenced 

1 

| w is an unsigned integer constant, less than or equal to 25 5, 

| specifying the number of characters of data 



The Z format code is used in conjunction with the transfer of hexade- 
cimal numbers. 

One storage location contains two hexadecimal digits. In read and 
write operations, padding and truncation are on the left. However, in a 
read operation, the padding character is a hexadecimal zero; in a write 
operation, it is a blank* 

If a :\ eight-byte internal field with the hexadecimal pattern 
* 01 23456789ABCDEF" is required, the external record could contain the 



characters 123456 7 8 9ABCDKF and would be read by a Z15 format code. The 
high-order zero is automatically provided as the padding character. 

L Format Code 

l . — . — ... — _ . _. — _. — _. ._„_________„_____„___„„—...._--. ...__- — - 1 

1 General Form | 

y . . __„„_. „___„„„____„__„_„„„„„{ 

aLw | 

I 

where a is optional and is an unsigned integer constant, less than | 

or equal to 2 55, used to denote the number of times the same | 

format code is repetitively referenced | 

I 
w is an unsigned integer constant less than or equal to 25 5, j 
specifying the number of characters of data j 



Logical variables may be read or written by means of the format code 
Lw. 

Cn input, the first T or F encountered in the next w characters of 
the input record causes a value of .TRUE. or -.FALSE., respectively, to 
be assigned to the corresponding logical variable. If the field w con- 
sists entirely of blanks, a value of . FALSE. is assumed. 

On output, a T or an F is inserted in the output record corresponding 
to the value of the logical variable in the I/O list. The single 
character is preceded by w - 1 blanks. 

A Format Code 



j General Form 

|. __ 



aAw J 

1 

where a is optional and is an unsigned integer constant, less than | 

or equal to 255, used to denote the number of times the same f 

format code is repetitively referenced | 

w is an unsigned integer constant less than or equal to 25 5, \ 

specifying the number of characters of data j 



The format code Aw is used to read or write alphameric data. i t w is 
equal to the number of characters corresponding to the length speoir'ic.i- 
tion of each item in the I/O list, w characters are read or written. 

On input , if w is less than the length specif i cat ion of each item in 
the I/O list, w characters are read and the remaining rightmost charac- 
ters in the item are replaced with blanks. If w is greater than the 
length specification, the number of characters equal to the difference 
between w and the length specification are skipped, beginning with the 
leftmost character, and the remaining characters are read. 

On output, if w is less than the length specification of the item in 
the I/O list, the printed line will consist of the leftmost w character:.-, 
of the item. If w is greater than the length specification, the pr.i nted 
line will consist of the characters right- justified in the field and 
will be preceded by blanks. Hieref ore , it is important to always alloc- 
ate enough storage area *:o handle the characters being written (see "Ihe 
Type Statements"). 
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Example 1 : Assume that the array ALPHA consists of one subscript param- 
eter ranging from 1 through 20, The following statements could be writ- 
ten to copy a record from one data set to another. 



10 FORMAT (20A4) 



READ (5,10) (ALPHA(I),I=1,20) 



WRITE (6,10) (ALPHACI> # I=l f 20) 



Explanation : The READ statement would cause 20 groups of characters to 
be read from the data set associated with data set reference number 5. 
Each group of four characters would be placed into the 20 storage loca- 
tions starting with ALPHA (1) and ending with ALPHAC20). The WRITE sta- 
tement would cause the 20 groups of four characters to be written on the 
data set associated with data set reference number 6. 

Example 2 : As another example, consider all the variable names in the 
list of the READ statement, below, to have been explicitly specified as 
REAL, and the array CONST to have been specified as having one subscript 
parameter ranging from 1 through 10. Then assume this input data is 
associated with data set reference number 5 

ABCDE. . .XYZ$1234567890b 

where . .. represents the alphabetic characters F through W, and b means 
a blank. These statements could be written 



10 FORMAT (27A1,10A1,A1) 
20 FORMAT C • x ' , 6 ( 7A1 , 5X) ) 



READ (5,10)A,B,C,D,E,F,G,H,I, 

1 J,K,L,M,N,0,P,Q,R r 

2 S f T,U f V r W,X,Y,Z f $, 

3 (CONST (IND),IND=1, 10), BLANK 



DO 50 INDEX = 1,5 



WRITE ( 6 , 20 ) G,R ,0 ,U ,P , BLANK , CONST ( I NDEX ) , 

1 B,L,0,C,K, BLANK, CONST (I NDEX), 

2 F,I,E,L,D,BLANK,CONST(INDEX) , 

3 G,R,C,U,P , BLANK, CONST (INDEX+5), 

4 B,L,0,C,K, BLANK, CONST (INDEX+5), 

5 F, I, E,L,D, BLANK, CONST( INDEX+5) 



50 CONTINUE 



50 



Explanation : The READ statement would cause the 37 alphameric charac- 
ters and the blank in the data set associated with data set reference 
number 5 to be placed in the storage locations specified by the variahle 
names in the READ list. Thus, the variables A through Z receive the 
values A through Z, respectively? the variahle $ receives the value $; 
the numbers 1 through 9, and 0, are placed in the 10 fields in storage 
starting with CONST(l) and ending with CONST(IO); and the variable BLANK 
receives a blank. The WRITE statement within the DO loop would cause 
the following heading to be printed. A subsequent WRITE statement 
within the DO loop could then be written to print the corresponding out- 
put data. 

Print Position 1 Print Position 67 

• ■ 

I I 

GROUP 1 BLOCK 1 FIELD 1 GROUP 6 BLOCK 6 FIELD 6 

(output data) - 

GROUP 2 BLOCK 2 FIELD 2 GROUP 7 BLOCK 7 FIELD 7 

(output data) - 



GROUP 5 BLOCK 5 FIELD 5 GROUP BLOCK FIELD 

(output data) - 

Literal Data in a Format Statement 

Literal data consists of a string of alphameric and special charac- 
ters written within the FORMAT statement and enclosed in apostrophes. 
The string of characters must be less than or equal to 255. For example: 

25 FORMAT C 1970 INVENTORY REPORT*) 

An apostrophe within the string is represented by two successive apo- 
strophes; for example, the characters DON'T are represented as: 

DON' § T 

The effect of the literal format code depends on whether it is used 
with an input or output statement. 

INPUT 

A number of characters, equal to the number of characters between the 
apostrophes, are read from the designated data set. These characters 
replace, in the FORMAT statement, the characters within the apostrophes. 
For example. 



FORMAT (• HEADINGS*) 



READ (3,5) 
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would cause the next, nine characters to be read from the data set asso- 
ciated with data set reference number 3; these characters would replace 
the blank and the eight, characters SHEADINGS in the FORMAT 
statement . 

C.HJTPUT 

All characters (including blanks) within the apostrophes are written 

• is part of the output data; thus. 



FORMAT (• THIS IS ALPHAMERIC DATA 1 ) 



WRITE (2,5) 



would «>iuse the following record to be written on the data set asso- 
.;idt*Hi with the data set reference number 2 

tfHiS IS ALPHAMERIC DATA 

*?i;i*rr I) indicates a blank. 
:\ hot w<M-_fJ\>de 

j iienr:ral Form | 

! w!« | 

t ! 

| wht-r. e w is an unsigned integer constant less than or equal to 255, 1 
| specifying the number of characters following H | 

Th^ H for ni at code is used in cod junction with the transfer of literal 

a , » * a . 

The tor mat code wH is followed in the FORMAT statement by w (w<255) 

\ , 1 1 a i : i e i s ; f o r example, 

L i » okMAT <31H THIS IS ALPHAMERIC INFORMATION) 

blanks ace significant and must be included as part of the count w. 
• s ' ih'> * of wH depends on whether it is used with input or output. 

1. < «n input , w characters are extracted from the input record and 
rtvi-'icc the w characters of the literal data in the FORMAT 

statement . 

?, :)v "lit p.it f the w characters following the format code are written 

r, j ,.irt of the output record. 



an unsigned integer constant less than or equal to 255, 
:i tying the number of blanks to be inserted on output or 
mmber of characters to be skipped on input 



When the wX Cw£25S) format code is used with a RED statement (i.e,, 
on input) , w characters are skipped before the next data item is read 
in. For example, if a card has six 10-column fields of integer quanti- 
ties, and the second quantity is not to be read, then 

5 FORMAT (110, 10X, 4110} 

may be used with the appropriate READ statement. 

When the wX format code is used with a WRITE statement (i.e., on out- 
put), w characters are left blank. Thus, the facility for spacing 
within a printed line is available. For example, 

10 FORMAT Cx 1 , 3CF6. 2 f 5X) ) 

may be used with an appropriate WRITE statement to print this line 

123.45bbbbb817«32bbbbb524-67bbbhb 

T Format Code 



j General Form j 

1 

1 where w is an unsigned integer constant less than or equal to 255, 

| specifying the position in a FORTRAN record where the 

| transfer of data is to begin 



Input and output may begin at any position by using the format code 
Tw (w<255). Only when tbe output is printed does the correspondence 
between w and the actual print position differ. In this case, because 
of the carriage control character, the print position corresponds to 
w~l, as in 

5 FORMAT CT4 0, s 1970 INVENTORY REPORT - T80 r f DECEMBER* Tl, 8 PART 
NO. 10095') 

The FORMAT statement above would result in this printed line 

Print Print Print 

Position 1 Position 39 Position 79 

n ■ m 

PART NO. 10095 196<* INVENTORY REPORT DECEMBER 

These statements 

5 FORMAT CT«*0, • HEADINGS') 



READ C3,5) 

would cause the first 39 characters of the input data to be skipped, and 
the next 9 characters would then replace the blank and the characters 
HEADINGS in the FORMAT statement. 

The T format code may be used in a FORMAT statement with any type of 
format code. For example, this statement is valid 

5 FORMAT CT100, F10.3, T50, E9. 3, Tl, * ANSWER IS 1 ) 
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Scale Factor - P 

The representation of the data, internally or externally , may be 
modified by the use of a scale factor followed by the letter P preceding 
a format code. 

The scale factor is defined for input and output as 

external quantity = internal quantity x 10**scale factor 

For input, when scale factors are used in a FORMAT statement, they 
have effect only on real data which does not contain an E or D decimal 
exponent. For example, if input data is in the form xx.xxxx and, it is 
to be used internally in the form . xxxxxx , the format code used to 
effect this change is 2PF7.4. 

INPUT 

As another example, consider this input data 

27bbb-93.2094bb-175.8041bbbb55.3647 
where b represents a blank. 
These statements 

5 FORMAT (I2,3F11.4) 



READ (6,5) K,A,B,C 

would cause these variables in the list to assume these values 

K : 27 B 2 -175. 80*41 

A : -93.209H C s 55. 3647 

These statements 

5 FORMAT (I2,1P3F11.4) 

READ (6,5) K,A,B,C 

would cause these variables in the list to assume these values 

K : 27 B : -17.580*4 

A : -9.3209 C : 5.5364 

These statements 

5 FORMAT (I2,-1P3F11. 4) 

READ (6,5) K,A,B,C 

would cause the variables in the list to assume these values 

K : 27 B : -1758. 041x 

A : -932.094x C : 553.647x 

where x represents an extraneous digit. 
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OUTPUT 

Assume that the variables K, A, B, and C have these values 

K j 27 B : -175.8041 

A t -93.2094 C s 55,3647 

these statements 

5 FORMAT (I2,1P3F11. 4) 

WRITE (4,5) K,A,B,C 

would cause the variables in the list to output these values 

K : 27 B j -1758. 041x 

A j -932.094x C : 553. 647x 

where x represents an extraneous digit.. 

These statements 

5 FORMAT (I2,-1P3F11.4) 

WRITE (4,5) K,A,B,C 

would cause the variables in the list to output these values 

K : 27 B : -17.5804 

A i -9.3209 C : 5.5364 

For output, when scale factors are used, they have effect only on real 
data. However, this real data may contain an E or D decimal exponent. 
A positive scale factor used with real data that contains an E or D 
decimal exponent increases the number and decreases the exponent. Thus, 
if the real data was in a form using an E decimal exponent, and the sta- 
tement FORMAT (IX, 12, 3E13.3) used with an appropriate WRITE statement 
resulted in this printed line 

27bbb-0.932Eb02bbb-0.175Eb03bbbb0.553Eb02 

then the statement FORMAT (1X,I2,1P3E13.3 ) used with the same WRITE sta- 
tement would result in this printed output 

27bbb-9.3 20Eb01bbb-1.758Eb0 2bbbb5. 536Eb01 

The scale factor is assumed to be zero if no other value has been 
given. However, once a value has been given, it will hold for all for- 
mat codes (i.e., those that correspond to real data) following the scale 
factor within the same FORMAT statement. This also applies to format 
codes enclosed within an additional pair of parentheses. Once the scale 
factor has been given, a subsequent scale factor of zero in the same 
FORMAT statement must be specified by OP. 

Carriage Control 

Vfrien records written under format control are prepared for printing, 
the following convention for carriage control applies: 

First Character Can lage Advance Before Printing 

Blank One line 

Two lines 

1 To first line of next page 
♦ No advance 
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Ihe first character of the output record may be used for carriage 
control and does not appear in the first print position of the print 
line. However, it appears in all other media as data. 

Carriage control can be specified in either of two forms of literal 
data* These statements would both cause two lines to be skipped before 

printing 

10 FORMAT PCT, 5CF7.3)) 
10 FORMAT (1H0, 5CF7.3)) 

ADDITIONAL INP UT/OUTP UT STA TEMENTS 

Hie statements END FILE, REWIND, and BACKSPACE are used to control 
the data sets, as described in the following text, 

END FILE Statement 

J General Form | 

i END FILE a | 

I I 

| where a is an unsigned integer constant or integer variable of | 
j length U that represents a data set reference number j 

L„. _____„„____. ______ . _ „ .__„_„____.„„. „_ .__ I 

The END FILE statement defines the end of the data set associated 
with a. A subsequent WRITE statement defines the beginning of a new 
data set, 

REWIND Statement 



r ._„__„_ ______ _, _. .__ _________ .___. . ___ j 

j General Form | 

K _-_„_-__„-_.„__„.__„„„__.._____„_. _ __ .j 

| REWIND a | 

1 I 

j where a is an unsigned integer constant or integer variable of | 

J length 4 that represents a data set reference number j 

i—_ ._______.. _____ ___ — — _ — . — __.. . — ____ — .__ — , — . — .. , — ___. „ . i 



The REWIND statement causes a subsequent READ or WRITE statement ref- 
erring to a to read data from or write data into the first data set 
associated with a. REWIND causes a logical rewinding to the beginning 
of the first data set associated with the specified data set reference 
number; it does not always cause a physical rewinding of the tape. If, 
however , the specified data set is the first on the tape, physical, as 
well as logical, rewinding occurs. 

BACKSPACE Statement 



r „__ — __„__ — ___. . — - — .___ . — . .__ — -__ . — ._. ._ , 1 

i General Form | 

j BACKSPACE a | 

i ! 

| where a is an unsigned integer constant or integer variable of | 

j length 4 that represents a data set reference number j 

L_____„___„__. ., . . . .__. , _ _. __J 

The BACKSPACE statement causes the data set associated with a to bac- 
kspace one record. If the data set associated with a is already at its 
beginning, execution of this statement has no effect, 
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SPECIFICATION STATEMENTS 



The specification statements provide the compiler with information 

about the nature of the data used in the source program. In addition, 
they supply the information required to allocate storage locations for 
this data. Specification statements describing data may appear anywhere 
in the source program, but must precede any statements which refer to 
that data. The specification statements are the type statements IMPLI- 
CIT, INTEGER, REAL, COMPLEX, and LOGICAL, and the DIMENSION, COMMON, and 
EQUIVALENCE statements. 

THE TYPE STATEMENTS 

There are two kinds of type statements i the IMPLICIT specification 
statement and the Explicit specification statements (INTEGER, REAL, COM- 
PLEX, and LOGICAL) . 

The IMPLICIT specification statement enables the user to 

1. Specify the type of a group of variables or arrays according to the 
initial character of their names, 

2. Specify the amount of storage to be allocated for each variable 
according to the associated type. 

The Explicit specification statements enable the user to 

1. Specify the type of a variable or array according to their particu- 
lar name. 

2. Specify the amount of storage to be allocated for each variable 
according to the associated type. 

3. Specify the dimensions of an array. 

U . Assign initial data values for variables and arrays. 
IMPLICIT Statement 

I General Form | 

IMPLICIT type »s C a lf a 2 , ...),... , type *s(a , ,a ? ,. . . ) 

where type represents one of the following: INTEGER, RKAL, COM- 
PLEX, or LOGICAL 

*s is optional and represents one of the permissible length 
specifications for its associated type 

i* i » f*2»-** represent single alphabetic characters each 

separated by commas, or a range of characters (in alphabetic 
sequence) denoted by the first and last characters of the 
range separated by a minus sign (e.g., (A-D) ) 



The IMPLICIT statement, if specified, should be the first statement 
in a main program, and the second statement in a FUNCTION, SUBROUTINE, 
or BLOCK DATA subprogram. 

The IMPLICIT type statement enables the user to declare the type of 
the variables appearing in his program (i.e., integer , real, complex, or 
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logical) by specifying that variables beginning with certain designated 
letters are of a certain type. Furthermore, the IMPLICIT statement 
allows the programmer to declare the number of locations to be allocated 
for each specified variable in a group. Hie types that a variable can 
assume, and the permissible lengths are 

Type Length Specification 

INTEGER 2 or U (standard length is 4) 

REAL 4 or 8 (standard length is 4) 

COMPLEX 8 or 16 (standard length is 8) 

LOGICAL 1 or 4 (standard length is 4) 

For each type there is a corresponding standard length specification. 
If this standard length specification (for its associated type) is 
desired, the *s may be omitted in the IMPLICIT statement. That is, the 
variables will assume the standard length specification. For each type 
there is also a corresponding optional length specification. If this 
optional length specification is desired, the *s must be included within 
the IMPLICIT statement. 

Example I s 

IMPLICIT REAL (A-H, 0-Z,$), INTEGER (I-N) 

Explanation : All variables beginning with the characters I through N 
are declared as INTEGER. Since no length specification was explicitly 
given (i.e., the *s was omitted), four storage locations (the standard 
length for INTEGER) are allocated for each variable. 

All other variables (those beginning with the characters A through H, 
O through Z, and $) are declared as REAL with four storage locations 
allocated for each. 

Note that the statement in Example 1 performs the same function of 
typing variables as the predefined convention (see "Type Declaration by 
the Predefined Specification"). 

Example 2 : 

IMPLICIT INTEGER*2(A-H), REAL*8(I-K) f LOGICAL (L,M ,N) 

Explanation : All variables beginning with the characters A through H 
are declared as integer, with two storage locations allocated for each. 
Ml variables beginning with the characters I through K are declared as 
real, with eight storage locations allocated for each. All variables 
beginning with the characters L f M, and N are declared as logical, with 
four locations allocated for each- 

Since the remaining letters of the alphabet (O through Z and $) were 
left undefined by the IMPLICIT statement, the predefined convention will 
take effect. Thus, all variables beginning with the characters 
through Z and $ are declared as real, each with a standard length of . 
four locations each. 

E xample 3 : 

IMPLICIT C0MPLEX*16(C-F) 

Explanation : All variables beginning with the characters C through F 
are declared as complex, each with eight storage locations reserved for 
the real part of the complex data and eight storage locations reserved 
for the imaginary part. The types of the variables beginning with the 
characters A, B, G through Z, and $ are determined by the predefined 
convention. 
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Explicit Specification Statements 

| General Form [ 

type *s a*s t (k, )/Xj/,b*s ? (k 2 )/x 2 /,.. . ,jz*s n Ck n )/x n / 

where Type is INTEGER , REAL, LOGICAL, or COMPLEX 

*s, *s l , *s 2 ,. . . ,*s n are optional? each £ represents one of 
the permissible length specifications for its associated 

type 

a,b,...,£ represent variable, array, or function names (see 
"SUBPROGRAMS") 

(k j) , (k 2 ) , . . , ,(k n ) are optionali each k is composed of 1 
through 7 unsigned integer constants, separated by commas, 
representing the maximum value of each subscript in the 
array? each k may be an unsigned integer variable only when 
it appears in a Type statement in a subprogram 

/x i/,/x 2 /,. . . #/x n / are optional and represent initial data 
values 



The Explicit specification statements declare the type (INTEGER, 
REAL, COMPLEX, or LOGICAL) of a variable or array by its name, rather 
than by its initial character. This differs from the other ways of spe- 
cifying the type of a variable or array (i.e., the predefined convention 
and the IMPLICIT statement). Also, the information necessary to alloc- 
ate storage for arrays (dimension information) may be included within 
the statement. However, if this information does not appear in an 
Explicit specification statement, it must appear in a DIMENSION or COM- 
MON statement (see "DIMENSION Statement" or "COMMON Statement"). 

Initial data values may be assigned to variables or arrays by use of 
/x n /, where x n is a constant or list of constants separated by commas. 
This set of constants may be in the form "r* constant", where r is an 
unsigned integer, called the repeat constant. 

No element may have more than one initial value given in the same 
program. A function name may not have an initial value assigned to it. 
An initially defined variable or a variable of an array may not be in 
blank common; in a labeled common block, they may be initially defined 
only in a BLOCK DATA subprogram. Initial data values may not be 
assigned to dummy segments appearing in a FUNCTION, SUBROUTINE, or ENTRY 
statement. 

In the same manner in which tbe IMPLICIT statement overrides the pre- 
defined convention, the Explicit specification statements override the 
IMPLICIT and predefined convention. If the length specification is 
omitted (i.e.,*s), the standard length per type is assumed. 

Example 1 : 

INTEGER*2 ITEM/76/, VALUE 

Explanation ; This statement declares that the variables ITEM and VALUE 
are of type integer, with two storage locations reserved for each. 
Also, the variable ITEM is initialized to the value 76. 

Example 2 : 

COMPLEX C,D/(2. 1,4.7)/, E*16 

Explanation : This statement declares that the variables C, D, and E are 
of type complex. Since no length specification was explicitly given for 
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C and D, the standard length is assumed. Thus, C and D have eight 
storage locations reserved for each (four for the real part, four for 
the imaginary part) and D is initialized to the value (2,1,4.7). In 
addition, 16 storage locations are reserved for the variable E. Thus f 
if a length specification is explicitly written, it overrides the 
assumed standard length. 

Example 3 ; 

REAL+8 ARRAY, HOLD, VALUE***, ITEM(5,5> 

Explanation s Ihis statement declares that the variables ARRAY, HOLD, 
VALUE, and the array named ITEM are of type real. In addition, it 
declares the size of the array ITEM, ARRAY and HOLD have eight storage 
locations reserved for each; VALUE has four locations reserved; and ITEM 
has 200 storage locations reserved (eight for each variable in the 
array). Note that when the length is associated with the type (e.g., 
REAL*8), the length applies to each variable in the statement unless 
explicitly overridden (as in the case of VALUE*4). 

Example 4 : 

REAL ACS, 5)/20*6.9E2, 5*1.0/, B(100 )/100*0 . 0/,TOAD*8 (5)/5* 0. 0/ 

Explanation : This statement declares the size of each array, A and B, 
and their type (real) . The array A has 100 storage locations reserved 
(four for each variable in the array); the array B has 400 storage loca- 
tions reserved (four for each variable). Also, the first 20 variables 
in the array A are initialized to the value 6.9E2 and the last five 
variables are initialized to the value 1. 0. All 100 variables in the 
array B are initialized to the value 0.0. The array TOAD has 40 storage 
locations reserved (eight for each variable). Also, each variable is 
initialized to the value 0.0. 

Example 5 : 

REAL A/Z1234CAF9/,B 

Explanation : This statement declares that the variables A and B are of 
type real, each with four storage locations reserved. Also, variable A 
is initialized to 1234CAF9 by using the hexadecimal constant. Note that 
the maximum number of digits allowed is dependent upon the length speci- 
fication of the variable being initialized. If the number of digits is 
greater than the maximum allowed, the leftmost hexadecimal digits are 
truncated; if less than the maximum, hexadecimal zeros are supplied on 
the left (see "Hexadecimal Constants"). 

ADDITIONAL SPECIFICATION STATEMENTS 
DIMENSION Statement 

| General Form | 

^_ 1 

DIMENSION a j (k t ) ,a 2 (k 2 ) , a 3 (k 3 ) ,. . . ,a n (k n ) 

where a,, a 2 , a 3 ,..., a n are array names 

JS i » iS?' k 3 ,.... ,k n are each composed of 1 through 7 unsigned 
integer constants, separated by commas, representing the 
maximum value of each subscript in the array; k { through k n 
may be integer variables of length 4 only when they appear 
in a DIMENSION statement within a subprogram. 
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Hie information necessary to allocate storage for arrays used in the 
source program may be provided by the DIMENSION statement . The follow- 
ing examples illustrate how this information may be declared. 

Examples ; 

DIMENSION AC10) f ARRAY (5,5,5,5,5), LISTC10,100) 
DIMENSION B(25,50) ,TABLE(25 ,25 ,25 ) 

Adjustable Dimensions 

The previous examples showed that the maximum value of each subscript 
in an array was specified by a numeric value. These numeric values 
(maximum value of each subscript) are known as the absolute dimensions 
of an array and may never be changed. However, if an array is used in a 
subprogram (see "Subprograms") and is not in Common, the size of this 
array does not have to be explicitly declared in the subprogram by a 
numeric value. That is, the Explicit specification statement, appearing 
in a subprogram, may contain integer variables that specify the size of 
the array. When the subprogram is called, these integer variables 
receive their values from the calling program. Thus, the dimensions 
(size) of a dummy array appearing in a subprogram are adjustable and may 
change each time the subprogram is called. 

The absolute dimensions of an array must be declared in a calling 
program. The adjustable dimensions of an array, appearing in a subpro- 
gram, should be less than or equal to the absolute dimensions of that 
array, as declared in the calling program. 

The following example illustrates the use of adjustable dimensions; 

Example : 

Calling Program Subprogram 



REAL*8 A(5,5) SUBROUTINE MAPMY (• . • , R,L,M, . . •) 



CALL MAPMYC* ,A,2,3,...) REAL*8« . • ,R (L,M) , . 



DO 100 1=1, L 



Explanation s Hie statement REAL*8 ACS, 5) appearing in the calling pro- 
gram declares the absolute dimensions of the array A. When the subrou- 
tine MAPMY is called, the dummy argument R assumes the array name A, and 
the dummy arguments L and M assume the values 2 and 3, respectively, 

Hie correspondence of the subscripted variables of the arrays A and R 
is shown in the following example. 




AC2,2) 



Specification Statements 61 



Thus, in the calling program, the subscripted variable A(l, 2) refers to 
the sixth subscripted variable in the array A. However, in the subpro- 
gram MAPMY the subscripted variable R(l,2) refers to the third sub- 
scripted variable in the array A, which is A(3 f l). This is so because 
the dijmens ions of the array R , as declared in the subprogram, are not 
the same as those in the calling program. 

If the absolute dimensions in the calling program were the same as 
the adjusted dimensions in the subprogram, the subscripted variables 
R<1,1) through R(5,5) in the subprogram would always refer to the same 
storage locations as specified by the subscripted variables AC 1,1) 
through A (5,5) in the calling program, respectively. 

The numbers 2 and 3, which became the adjusted dimension of the dummy 
array R, could also have been variables in the argument list of the cal- 
ling program. For example, assume that this statement was in the cal- 
ling program: 



CALL MAPMY C...,A,I,J,, 



) 



Then, as long as the values of I and J were previously determined, the 
arguments may be variables. Also, the variable dimension size may be 
passed through more than one level of subprograms; for example, within 
the subprogram MAPMY could have been a call statement to another subpro- 
gram in which dimension information about A could have been passed. 

If an array has a variable dimension, that array name must be a dummy 
variable (i,e., must appear in a FUNCTION, SUBROUTINE, or ENTRY state- 
ment) . The variable dimension itself can be a dummy variable or can 
appear in a COMMON statement. 



COMMON Statement 



j General Form 
| 

COMMON 



^ 



/r/a (k,),b(k 2 ), /£/c(k 3 ) ,d(k ),... 

where a,b, . . . ,c,d. . . are variable or array names 

k l# k 2# ,.,k 3 ,k ,.- are optional and are each composed of one 
through seven unsigned integer constants, separated by com- 
mas, representing the maximum value of each subscript in the 
array 

/r/..- represent optional common block names consisting of 
one through six alphameric characters, the first of which is 
alphabetic. These names must always be embedded in slashes 

. . . . . . j 



Variables or arrays that appear in a calling program or a subprogram 
may be made to share the same storage locations with variables or arrays 
in other subprograms by use of the COMMON statement. For example, if 
one program contains the statement 

COMMON TABLE 
and a second program contains the statement 

COMMON LIST 
the variable names TABLE and LIST refer to the same storage locations. 
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If the main program contains the statements! 



REAL A,B # C 

COMMON A,B f C 

and a subprogram contains the statements: 

REAL X,Y,Z 

COMMON X,Y,Z 

A shares the same storage location as X; B shares the same storage loca- 
tion as Y; and C shares the same storage location as Z. 

Consider the following examples: 

Example 1 : 



railing Program 



Subprogram 

SUBROUTINE MAPMY (...) 



COMMON A, B, C, RC10O) 

REAL A,B,C 

INTEGER R 



COMMON X, Y, Z, S(100) 
REAL X,Y,Z 
INTEGER S 



CALL MAPMY (, 



,) 



explanation : In the calling program, the statement COMMON A, B, C,R(1Q0) 
would cause 412 storage locations (four locations per variable) to be 
reserved in this manner 



Beginning 
of common 
area 



l___. 



_ — __ . j . — ._. . T f 

A | B 1 C | Layout of 

4 locations | ** locations | 4 locations | storage 
I I I 

_ ._. — _ + „ — .._._. + ~ .. 4 

R(l) I • . . 1 RC100) | 

4 locations | j ** locations j 

.. x . — . __. j . J 



The statement COMMON X, Y, Z, S(100) would then cause the variables 
X, Y # Z, and SCI) . . .S (100) to share the same storage space as A, B, C, 
and R(l) . . .R(100) . 

The example above shows that COMMON statements may be used to func- 
tion as a medium to implicitly transmit data from the calling program to 
the subprogram. That is, values for X, Y, Z, and S (1) . . . S(IOO) , because 
they occupy the same storage locations as A # B, C, and R(l ) . . .R (1 00) r do 
not have to be transmitted in the argument list of a CALL statement- 
Arguments passed through COMMON must follow the same rules of presenta- 
tion with regard to type, length, etc., as arguments passed in a list 
(see "SUBPROGRAMS") . 

Example 2 : Assume COMMON is defined in a main program and three subpro- 
grams as 

Main program — COMMON A,B,C 

Subprogram 1 — COMMON D,E,F 

Subprogram 2 -- COMMON Q,R,S,T,U 

Subprogram 3 -- COMMON V,W,X,Y,Z 

Also, assume the length specifications of these variables are so defined 
that the common area is shared as follows 
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I c 

I 2 locations 

1 r" 

( 2 locations 

i c I 

•if-. ' / L-.'-'rs?- i ons j 2 locations 

! Y I Z 

.•i;-; f ' ocii i ons j 2 locations 

") . , '• •::.,."'. i" validly referred to 
..,/!> i'\*r.o V,W,X,Y,Z. In addi- 

► C ( b , 'j t ,tri.* t 2;, It is also poss- 

i i -fj j a 'id R in Subprogram 2. 

v*' \'\->> : *c-'J in certain cases nay 
■. •<*..:! c ! ,, and y and H are 
if . »fid i:ua^ irtar y parts of D, 

■ • •.- I <■■; em:*- by the programmer must 



■?.' n ''in : : t. or age area ( coin- 
<i;-ii. Tnat is, no par- 
1 1 le va i: i a b 1 e s th a t 

locations relative to 
>r r variables and arrays 
• th,v jt i .separate areas 
j ! •> u q I i \: i x a lpha me r i c 
: : ■ c . r i c c.y s t h at h av e 



lamed ) common are pre- 
Imi example, the 

common area HOLD by 



,; . •; . f r 'in r.iished from laoeled 
c • 'i : by t wn consecuti ve 
1 niina ui" the common sta- 



,i«ie, iii 



. • i *i f " i :d in blank c ommon i n 

"< : ,• 1 a c ed. in the c o mmo n area 



Kti • :*; ;. is COMMON statements are 



statement 



Example 3 : Assume that A, B, C, K, X, and Y each occupy four locations 
of storage, H and G each occupy eight locations, and D and E each occupy 
two locations. 



Calling Program 



Subprogram 
SUBROUTINE MAPMYU..) 



COMMON H, A /R/ X, D // B 



COMMON G, Y f C /R/ K, E 



CALL MAPMY ( • • • ) 



Explanation : In the calling program, the statement COMMON H,A/R/X,D//B 
causes 16 locations (four locations each for A and B, and eight for H) 
to be reserved in blank common in this order 



Beginning ( 

of blank | 

common | 

h 

I 



H 

8 locations 



A B 

4 locations 4 locations 



1 



continuation of blank common 






and also causes six locations (four for X and two for D) to be reserved 
in this labeled common area R in this order 



Beginning j 
of labeled! 
common R | 

h 
I 



X D 

U locations 2 locations 



-H 



continuation of labeled common 



l_ 



The statement COMMON G,Y,C/R/K,E appearing in the subprogram MAPMY 
would cause the variables G, Y, and C to share the same storage space 
(in blank common) as H, A, and B f respectively. It would also cause the 
variables K and E to share the same storage space (in labeled common 
area R) as X and D, respectively. The length of a COMMON area may be in- 
creased by using an EQUIVALENCE statement (see "EQUIVALENCE Statements") 

Programming Consideration s 

Variables in a COMMON block may be in any order. However, consider- 
able object -time efficiency is lost unless the programmer ensures that 
all the variables have proper boundary alignment. 

Proper alignment is achieved either by arranging the variables in a 
fixed descending order according to length, or by constructing the block 
so that dummy variables force proper alignment. If the fixed order is 
used, the variables must appear in this order 

length of 16 (complex) 

length of 8 (complex or real) 

length of 4 (real or integer or logical) 

length of 2 (integer) 

length of 1 (logical) 
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Type of 
Variable 


Length 
Speci f ica tion 


Logical 


1 

4 


Integer 


2 


Real 


8 


Complex 


8 
16 



If the fixed order is not used, proper alignment can be ensured by 
constructing the block so that the displacement of each variable can be 
evenly divided by the reference number associated with the variable. 
(Displacement is the number of storage locations from the beginning of 
the block to the first storage location of the variable.) The following 
list shows the reference number for each type of variable. 

Reference 

Number 
1 

2 
I* 

8 

8 
8 

The first variable in every COMMON block is positioned as if its 
length specification were eight. Therefore, a variable of any length 
may be the first assigned within a block. To obtain the proper align- 
ment for other variables in the same block, it may be necessary to add a 
dummy variable to the block. For example, the variables A, I, and CMPLX 
are REAL*4, INTEGER* 4, and C0MPLEX*8, respectively, and form a COMMON 
block that is defined as 

COMMON A, I, CMPLX 

Then, the displacement of these variables within the block is illus- 
trated by 

A I CMPLX 

4 storage ** storage 8 storage 

locations locations locations 

displacement displacement displacement 

storage 4 storage 8 storage 

locations locations locations 

The displacements of I and CMPLX are evenly divisible by their 
reference numbers. However, if I were an integer with a length specifi- 
cation of 2, then CMPLX is not properly aligned (its displacement of 6 
is not evenly divisible by its reference number of 8). In this case, 
proper alignment is ensured by inserting a dummy variable with a length 
specification of 2 either between A and I or between I and CMPLX. 

EQUIVALENCE Statement 



| General Form | 

EQUIVALENCE (a, b, c, ), (d, e, f, ) 

where a, b, c, d, e, f,... are variables that may be subscripted. 
The subscripts may have two forms: If the variable is sing- 
ly subscripted it refers to the position of the variable in 
the array (i.e., first variable, 25th variable, etc); if 
the variable is multi subscripted, it refers to the position 
in the array in the same manner as the position is referred 
to in an arithmetic statement 
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Hie EQUIVALENCE statement provides the option for controlling the 
allocation of data storage within a single program or subprogram. It is 
analogous to the option of using the COMMON statement to control the 
allocation of data storage among several programs. When the logic of 
the program permits, the number of storage locations used can be reduced 
by causing locations to be shared by two or more variables of the same 
or differing types and lengths. The EQUIVALENCE statement cannot be 
used to obtain mathematical equality of two variables. 

Example 1 : 

DIMENSION B(5) f C(10, 10) r D(5 , 10, 15) 
EQUIVALENCE (A, B(l), C(5,3)), (D(5,10,2), E) 

Explanation : This EQUIVALENCE statement indicates that the variables 
A, BCD, and CCS, 3) are assigned to the same storage locations; also that 
DC 5, 10, 2) and E are assigned to the same storage locations. In this 
case, the subscripted variables refer to the position in an array in the 
same manner as the position is referred to in an arithmetic statement. 
Note: Variables or arrays that are not mentioned in an EQUIVALENCE sta- 
tement are assigned to unique storage locations. The EQUIVALENCE state- 
ment must not contradict itself or any previously established equiva- 
lences. For example, the further equivalence specification of BC 2) with 
any other element of the array C, other than C(6,3), is invalid. 

Example 2 : 

DIMENSION B(5), CC10, 10), DCS, 10, 1S>< 
EQUIVALENCE (A, BCD, CC2S)), (DC100) , E) 

Explanation : This EQUIVALENCE statement indicates that the variable A, 
the first variable in the array B, namely BCD, and the 25th variable in 
the array C, namely C(5,3), are to be assigned the same storage loca- 
tions. Also, it also specifies that D(IOO), i.e., DC5,10,2), and E are 
to share the same storage locations. Notes The effects of the EQUIVA- 
LENCE statements in examples 1 and 2 are the same. 

Variables that are brought into COMMON through EQUIVALENCE statements 
may increase the size of the block, as indicated by these statements: 

COMMON A, B, C 
DIMENSION DC3) 
EQUIVALENCE (B,D(1)> 

This would cause a common area to be established containing the 
variables A, B, and C. The EQUIVALENCE statement would then cause the 
variable D(l) to share the storage location with B, D(2) would share 
with C, and D(3) would extend the size of the common area, in this 
manner 

A (lowest location of the common area) 

B, D(l) 

C, DC2) 

D(3) (highest location of the common area) 

Since arrays must be stored in consecutive forward locations, a vari- 
able may not be made equivalent to another variable of an array in such 
a way as to cause the array to extend before the beginning of the common 
area. For example, this EQUIVALENCE statement is invalid 

COMMON A, B, C 
DIMENSION D(3) 
EQUIVALENCE (B, D(3)> 
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because it would force D(l> to precede A, as follows: 

DC1> 

A, DC 2) (lowest location of the common area) 

B, DC3) 

C (highest location of the common area) 

Programming Considerations 

Two variables in one COMMON block or in two different COMMON blocks 
may not be made equivalent. Variables in an equivalence group may be in 
any order . However, considerable object-time efficiency is lost unless 
the programmer ensures that all the variables have proper boundary 
alignment. 

Proper alignment is achieved either by arranging the variables in a 
fixed descending order according to length, or by constructing the group 
so that dummy variables force proper alignment. If the fixed order is 
used, the variables must appear in this order: 

length of 16 (complex) 

length of 8 (complex or real) 

length of 4 (real or integer or logical) 

length of 2 (integer) 

length of 1 (logical) 

If the fixed order is not used, proper alignment can be ensured by 
constructing the group so that the displacement of each variable in the 
group can be evenly divided by the reference number associated with the 
variable- (Displacement is the number of storage locations from the 
beginning of the group to the first storage location of the variable.) 
(The reference numbers for each type of variable are given under "COMMON 
Statement.") The first variable in each group is positioned as if its 
length specification were eight. 

For example, the variables A, I, and CMPLX are REAL* 4, INTEGER***, and 
COMPLEX*8, respectively, and are defined as 

DIMENSION A(10), 1(16), CMPLX(5) 
EQUIVALENCE (A(l), 1(7), CMPLX(D) 

Hien, the displacement of these variables within the group is illus- 
trated by 



1(1) 



64 storage locations 
40 storage locations 
40 storage locations 



I (16) 
A (10) 
CMPLX (5) 



AC1) 

CMPLX (1) 

displacement displacement 
storage 24 storage 

locations locations 

The displacements of A and CMPLX are evenly divisible by their reference 
numbers. However, if the EQUIVALENCE statement were written as 

EQUIVALENCE (A(l), 1(6), CMPLX (1)) 

then CMPLX is not properly aligned (its displacement of 20 is not evenly 
divisible by its reference number, 8). 
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NAMING SUBPROGRAMS 



A subprogram name coaf.isrr 
ters, t he f ir s t o f wh ic h mi s v 
contain special characters i~ 
may be typed Cas variables are) 



Type declaratio: \ 
accomplished in on: ■ > 
by the .IMPLICIT statement . . -; 
ments . Thus, the san*? r a »*.":> 1 
app ly t o Sta t ement Fu not i - vc "» . 



Type declara tion 01 _ FyNc/i 
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IMPLICIT statement, or i.y .: 
cification of the FUNCTION 



3 . Type declaration of 



u bi 



1 -a? 



TINE subprogram cannot. 

returned to the call) a: r, *i 
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A function is a s t.a te me n t 
variables. To use a fuiiov 101 



1. Define the function (i.e., specify what calculations are to be 
performed) , 

2. Refer to the function by name* where required in the program. 
Function Definition 

The three steps in the definition of a function are 

1. The function must be assigned a unique name by which it may be 
called (see "Naming Subprograms" K 

2. The arguments of the function must be stated, 

3. The procedure for evaluating the function must be stated. 

Items 2 and 3 are discussed in detail in the sections dealing with 
the specific subprogram ("Statement Functions," "FUNCTION Subprograms," 
etc, ) . 

Function Reference 

The name of a function, appearing in any FORTRAN arithmetic expres- 
sion, refers to the function. Thus, the appearance of a function, with 

its arguments in parentheses, causes the computations to be performed as 
indicated by the function definition* The resulting quantity replaces 
the function reference in the expression and assumes the type of the 
function. The type and length of the name used for the reference must 
agree with the type and length of the name used in the definition. 

STATEMENT FUNCTIONS 

Statement functions are defined by a single arithmetic or logical 
assignment statement within the program in which they appear. For 
example, 

FUNC(A,B) = 3.*A+B**2.+X+Y+Z 

defines the statement function FUNC, where FUNC is the function name and 
A and B are the function arguments. 

The expression on the right defines the computations that are to be 
performed when the function is used in an arithmetic statement. This 
function might be used in this way 

C = FUNC(D,E) 

which is equivalent to 

C = 3,*D+E**2.+X+Y+Z 

Nate the correspondence between A and B in the function definition 
statement and D and E in the arithmetic statement. The quantities A and 
B, enclosed in parentheses following the function name, are the argu- 
ments of the function. They are dummy variables for which the quanti- 
ties D and E, respectively are substituted when the function is used in 
an arithmetic statement. 
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I - — ■ ■ • ■ 1 

I General Form | 

name (a,b, ...,n) = expression 

where name is any subprogram name (see "Naming Subprograms" ) . 

a,b,...,n are distinct (within the same statement) nonsub- 
scripted variables 

expression is any arithmetic or logical expression that does 
not contain subscripted variables; any statement functions 
appearing in this expression must be defined previously 



The actual arguments must correspond in order, number, and type to 
the dummy arguments. At least one argument must be used, A maximum of 
15 variables appearing in the expression may be used as arguments of the 
fun ct ion . 

Note: All Statement Function definitions to be used in a program must 
precede the first executable statement of the program. 

Examples : Valid statement function definitions: 

SUM(A,B,C,D) « A+B+C+D 
FUNC(Z) - A+X*Y*Z 
AVG(A,B,C,D) = (A+B+C+D)/** 
R0OT(A,B,C) = SQRT(A**2+B**2+C**2) 
VALID(A,B) = .NOT. A. OR. B 

Note : The same dummy arguments may be used in more than one Statement 
Function definition and, as variables, outside Statement Function 
definitions. 

Invalid statement function definitions 

SUBPRG (3, J f R)=3*I + J**3 (arguments must be variables) 

SOMEF(ACI) f B)=A(I)/B+3. (arguments must be nonsub- 

scripted) 
SUBPROGRAM (A, B) =A**2+B**2 (function name exceeds limit 

of six characters) 
3FUNC (D)=3.14*E (function name must begin with 

an alphabetic character) 
ASF(A)=A+B(I) (subscripted variable in the 

expression) 

Valid statement function references 

NET = GROS - SUM (TAX, FICA, HOSP, MISC) 

ANS = FUNC (RESULT) 

GRADE = AVG(LAB f LECTUR , SUMCTEST1, TEST2, TEST3, TESTU ) , FACTOR) 

Invalid statement function references 

WRONG = SUM (TAX, FICA) (number of arguments 

does not agree with 

above definition) 
MIX = FUNC (I) (mode of argument 

does not agree with 

above definition) 
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FUNCTION SUBPROGRAMS 

The FUNCTION subprogram is a FORTRAN subprogram consisting of any 
number of statements. It is an independently written program that is 

executed wherever its name appears in another program* 



General £orm 



FUNCTION name Ca 1# a 2 ,a 3 , ,a n ) 



f 



RETURN 



END 



where name is subprogram name (see "Naming Subprograms 1 *) 

aj*a 2f a3 r ««« # a n are n on subscripted variable, array, or dummy 
names of SUBROUTINE or other FUNCTION subprograms. There 
must be at least one argument in the argument list. (Argu- 
ments in a FUNCTION or SUBROUTINE subprogram may be enclosed 

in slashes within the commas. This form is equivalent to 
the normal format without the slashes, ) 

Since the FUNCTION is a separate subprogram, the variables and state- 
ment numbers within it do not relate to any other program. 

The FUNCTION subprogram may contain any FORTRAN statement except a 
SUBROUTINE statement, another FUNCTION statement, or BLOCK DATA 

statement . 



The arguments of the FUNCTION subprogram (i.e., a 1 ,a 2l a 3 ,...,a n ) may 
be considered as dummy variable names* These are replaced at the time 

of execution by the actual arguments supplied in the function reference 
in the calling program. The actual arguments may be: any type of con- 
stant, any type of subscripted or nonsubscrip ted variable, an array 
name, an arithmetic or logical expression, or the name of another sub- 
program. The actual arguments must correspond in number, order, type, 
and length to tbe dummy arguments. The array size must also be the 
same, except when adjustable dimensions are used. If the actual argu- 
ment corresponds to a dummy argument that is defined or redefined in the 
subprogram, the argument must be a variable name, subscripted variable 
name, or array name* All arguments in a subprogram refer to the storage 
area assigned to the arguments by the calling program. 

The relationship between variable names used as arguments in the cal- 
ling program and the dummy variables used as arguments in the FUNCTION 

subprogram is illustrated by 



Examp le 1 : 

Calli ng Program 

A = SOMEFCB,C) 



FUNCTION Subprogram 

FUNCTION SOMEFCX,Y) 
SOMEF = X/Y 
RETURN 
END 
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Explanatio n; The value of the variable B of the calling program is used 
in the subprogram as the value of the dummy variable X; the value of C 
is used in place of the dummy variable Y. Thus if B = 10.0 and C = 5.0, 
then A = B/C, which is equal to 2.0. 

The name of the function must be assigned a value at least once in 
the subprogram as the argument of a CALL statement, as a DO variable, as 
the variable name on the left side of an arithmetic statement, or in an 
input list (READ statement) within the subprogram. 

Example 2 : 

Calling Program FUNCTION Subprogram 

FUNCTION CALC (A,B,J> 



ANS = ROOTl*CALC(X,Y,I> 



I = J*2 



CALC = A**I/B 



RETURN 
END 

Explanation : The values of X, Y, and I are used in the FUNCTION subpro- 
gram as the values of A, B, and J, respectively. The value of CALC is 
computed, and this value is returned to the calling program, where the 
value of ANS is computed. Hie variable I in the argument list of CALC 
in the calling program is not the same as the variable I appearing in 
the subprogram. 

When a dummy argument is an array name, an appropriate DIMENSION or 
Explicit specification statement must appear in the FUNCTION subprogram. 
None of the dummy arguments may appear in an EQUIVALENCE or COMMON sta- 
tement within the subprogram nor may they be given initial data values. 

Type Specification of the FUNCTION Subprogram 

In addition to the three ways of declaring the type of a FUNCTION 
name (i.e., predefined convention, IMPLICIT statement, Explicit specifi- 
cation statement), there is the option of explicitly specifying the type 
of a FUNCTION name within the FUNCTION statement. 



j General Form j 

type FUNCTION name *s (a x ,a 2 ,a 3 ,. . . ,a n ) 

where type is INTEGER., REAL, COMPLEX, or LOGICAL 

name is the name of the FUNCTION subprogram 

*s is optional and represents one of the permissible length 
specifications for its associated type 

a ir d 2f a 3f ... f a n are nonsubscripted variable, array, or dummy 
names of SUBROUTINE or other FUNCTION subprograms. (There 
must be at least one argument in the argument list) 
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Example 1 : 

REAL FUNCTION SOMEF (A,B) 

SOMEF = A**2 ♦ B**2 

RETURN 
END 

Example 2 : 

INTEGER FUNCTION CALC*2 (X,Y,Z> 

CALC = X+Y+Z**2 



RETURN 
END 



Explanation : The FUNCTION subprograms SOMEF and CALC in examples 1 and 
2 are declared as type REAL (length 4 ) and INTEGER (length 2), 
respectively. 

RETURN and END Statements in a JUNCTION Subprogram 

All FUNCTION subprograms must contain both an END and at least one 
RETURN statement. The END statement specifies, for the compiler, the 
end of the subprogram? the RETURN statement signifies a logical conclu- 
sion of the computation and returns any computed value and control to 
the calling program. More than one RETURN statement may be used in a 
FORTRAN subprogram. 



Example : 



10 



FUNCTION DAV (D,E,F) 
IF (D-E) 10, 20, 30 
A = D+2.0*E 



5 A = 



20 DAV 



F+2.0*E 



= A*B**2 



RETURN 
30 DAV = B**2 



RETURN 
EI® 
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Multiple Entry into a FUNCTION Subprogram 

The standard entry into a FUNCTION subprogram is made by a function 
reference in an arithmetic expression, when the function reference uses 
the name defined in the FUNCTION statement. Entry is made at the first 
executable statement following the FUNCTION statement. 

It is also possible to enter a FUNCTION subprogram by a function 
reference to a name defined in an ENTRY statement in the FUNCTION sub- 
program. Entry is made at the first executable statement following the 
ENTRY statement. The name given in the FUNCTION statement is used to 
return the value of the function to the point of reference, rather than 
the name of the ENTRY statement. 



SUBROUTINE SUBPROGRAMS 

The SUBROUTINE subprogram is similar to the FUNCTION subprogram in 
many respects: the rules for naming FUNCTION and SUBROUTINE subprograms 
are the same, they both require an END statement, and they both contain 
the same sort of dummy arguments. Like the FUNCTION subprogram, the 
SUBROUTINE subprogram is a set of commonly used computations, but it 
does not need to return any results to the calling program, as does the 
FUNCTION subprogram. 

The CALL statement (discussed later in this section) is used in a 
main program or another subprogram to invoke a SUBROUTINE subprogram. 

Since the SUBROUTINE is a separate subprogram, the variables and sta- 
tement numbers within it do not relate to any other program. 

The SUBROUTINE statement must be the first statement in the subpro- 
gram. The SUBROUTINE subprogram may contain any FORTRAN statement 
except a FUNCTION statement, another SUBROUTINE statement, or a BLOCK 
DATA statement. If an IMPLICIT statement is used in a SUBROUTINE sub- 
program, it must immediately follow the SUBROUTINE statement. 



j General Form | 

h __ _ < 

SIBROUTINE name (a , ,a 2 , a 3 ,. . . ,a n ) 



RETURN 

END 

where name is the subprogram name (see •'Naming Subprograms") 

a 1 #a 2 ,a 3 , . . . ,a n are arguments. (There need not be any.) 
Each argument used must be a nonsubscripted variable or 
array name, the dummy name of another SUBROUTINE or FUNCTION 
subprogram, or of the form * where the character "*" denotes 
a return point specified by a statement number in the cal- 
ling program 

The SUBROUTINE subprogram may use one or more of its arguments to 
return values to the calling program. Any arguments so used must appear 
on the left side of an arithmetic statement or in an input list within 
the subprogram, as arguments of a CALL statement, as DO variables, or as 
arguments in a function reference. The SUBROUTINE name must not appear 
in any other statement in the SUBROUTINE subprogram. 
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The arguments (a if a 2 , a 3 ,...,a n ) may be considered as dummy variable 
names that are replaced at the time of execution by the actual arguments 
supplied in the CALL statement. The actual arguments must correspond in 
number, order, type, and length to the dummy arguments. The array size 
must also be the same except when adjustable dimensions are used. Dummy 
arguments may not appear in an EQUIVALENCE or COMMON statement within 
the subprogram nor may they be given initial data values. 

Example : The relationship between variable names used as arguments in 
the calling program and the dummy variables used as arguments in the 
SUBROUTINE subprogram is illustrated in this example. The object of the 
subprogram is to copy one array directly into another. 

Main Program SUBROUTINE Subprogram 

DIMENSION XC100) f Y(100) 

SUBROUTINE COPY(A,B,N) 

DIMENSION A (100) ,B (100) 

DO 10 I = 1, N 
CALL COPY (X,Y,K) 10 BCD = A CD 

RETURN 

END 

CALL Statement 

The CALL statement is used only to call a subroutine subprogram. 



| General Form | 

c Ai»l* nam€ ( a _ ]f a 2# a 3# ...,a n ) 

where name is the subroutine's subprogram name, or a name defined 
in an ENTRY statement in the SUBROUTINE subprogram 

a t ,a 2 ,a 3 , . . . ,a n are the actual arguments that are being sup- 
plied to the subroutine subprogram; each may be of the form 
€n where n is a statement number (see "RETURN Statements in 
a SUBROUTINE Subprogram") 



Hie standard entry into a SUBROUTINE subprogram is made by a CALL 
statement that refers to that subroutine * s subprogram name. Entry is 
made at the first executable statement following the SUBROUTINE state- 
ment. Also, it is possible to enter a SUBROUTINE subprogram by a CALL 
statement that refers to a name defined in an ENTRY statement in the 
SUBROUTINE subprogram. The EOTRY statement is described below. 

Examples : 

CALL OUT 

CALL MATMPY (X , 5, 40 ,Y,7 , 2) 
CALL QDRTIC ( X, Y, Z,ROOTl , R00T2 ) 
CALL SUB1 (X+Y+S^ABDF" ,SINE) 

The CALL statement transfers control to the subroutine subprogram and 
replaces the dummy variables with the value of the actual arguments that 
appear in the CALL statement. The arguments in a CALL statement may be: 
any type of constant, any type of subscripted or non sub scripted vari- 
able, an arithmetic expression, the name of a subprogram, or a statement 
number (see "RETURN Statements in a SUBROUTINE Subprogram"). 

The arguments in a CALL statement must agree in number, order, and 
type with the corresponding arguments in the subroutine subprogram. The 
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array sizes must also be the same in the subroutine and the calling pro- 
grams, except when adjustable dimensions are used (see "Adjustable 
Dimensions"). If an actual argument corresponds to a dummy argument 
that is defined or redefined in the referenced subprogram, the actual 
argument must be a variable name, subscripted variable name, or array 
name. All arguments in a subprogram refer to the storage area assigned 
to the arguments by the calling program. 



RETURN Statement in a SUBROUTINE Subprogram 



| General Form 
RETURN 
RETURN i 



— 1 
I 



where i is an integer constant or variable of length 4 whose 

value, say n, denotes the nth statement number in the argu- 
ment list of a SUBROUTINE statement 






The normal sequence of execution following the RETURN statement of 
a SUBROUTINE subprogram is to the next statement following the CALL in 
the calling program- It is also possible to return to any numbered 
statement in the calling program by using a return of the type where i 
is an integer constant or variable. Returns of the type RETURN may be 
made in either a SUBROUTINE or FUNCTION subprogram (see, "RETURN and 
END Statements in a FUNCTION Subprogram"). Returns of the type RETURN 
i^ may only be made in a SUBROUTINE subprogram. In a main program, a 
RETURN statement performs! the same function as a STOP statement. 



Example : 

Calling Program 



Subprogram 

SUBROUTINE SUB (X,Y,Z,*,*> 



10 CALL SUB <A,B,C,fc30,£40> 
20 Y = A + B 



30 Y = 



40 Y = 



A * C 



B + C 



100 IF CR) 200,300,400 
200 RETURN 
3 00 RETURN 1 
400 RETURN 2 
END 



END 

Explanation : Execution of statement 10 in the calling program causes 
entry into subprogram SUB. When statement 10 is executed, the return 
to the calling program will be to statement 2 0, 30, or 40, if R is 
less than, egual to, or greater than zero, respectively. 

A CALL statement that uses a RETURN i form may be best understood 
by comparing it to a CALL and computed GO TO statement in sequence. 
For example, 

CALL SUB (P,t20,Q,€35,R,t22) 
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is equivalent to: 

CALL SUB (P,Q,R,I) 
GO TO (20, 35, 22), I 

where the index I is assigned a value of 1 , 2 , or 3 in the called 

subprogram. 

ENTRY Statement 

The standard (normal) entry into a SUBROUTINE subprogram from the 
calling program is made by a CALL statement that references the sub- 
program name. The standard entry into a FUNCTION subprogram is made 
by a function reference in an arithmetic expression. Entry is made at 
the first executable statement following the SUBROUTINE or FUNCTION 
statement. 

A subprogram may also be entered (either SUBROUTINE or FUNCTION) by 
a CALL statement or a function reference that references an EOTRY sta- 
tement in the subprogram. Entry is made at the first executable sta- 
tement following the ENTRY statement. 



-i 



| General Form 

ENTRY name (a x ,a 2 ,a 3 f . . . ,a n ) 

where name is the name of an entry point containing from one to 

six alphabetic and/or numeric characters, the first of which 
is alphabetic 

a } ,a 2 ,a 3 , . . . ,a n are the dummy arguments corresponding to an 
actual argument in a CALL statement or in a function 
ref errence 

ENTRY statements do not affect control sequencing during normal 
execution of a subprogram. The order, type, and number of arguments 
need not agree between the SUBROUTINE or FUNCTION statement and the 
ENTRY statements, nor must the EOTRY statements agree among themselves 
in these respects. Each CALL or function reference, however, must agree 
in order, type, and number with the SUBROUTINE, FUNCTION, or ENTRY sta- 
tement that it references. Entry may not be made into the range of a 
DO; further, a subprogram may not reference itself directly or through 
any of its entry points. This statement is regarded as nonexecutable 
within its subprogram, . If it appears in a function subprogram the 
name given in the FUNCTION statement is still used to return the value 
of the function to the point of reference, rather than the name of the 
ENTRY statement. 

Example 1 : 

Calling Program Subprogram 

SUBROUTINE SUB1 (U, V,W,X, Y, Z) 

■ * 

1 CALL SUB1 (A,B,C,D,E,F) 

U = V 

2 CALL SUB2 (G,H,P) 

ENTRY SUB2 (T,U,V) 
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3 CALL SUB3 



ENTRY SUB3 



END 



Explanation : The execution of statement 1 causes entry into SUBl, 
starting with the first executable statement of the subroutine. Execu- 
tion of statements 2 and 3 also causes entry into the called program, 
starting with the first executable statement following the ENTRY SUB2 (T, 
U,V) and ENTRY SUB3 state me nts, respectively. 

Entry into a subprogram initializes all references in the entire sub- 
program to items in the argument list. Return from a subprogram is made 
by way of the entry point referenced. ENTRY statements may only appear 
in JUNCTION or SUBROUTINE subprograms. The dummy arguments in a subpro- 
gram may appear in any statement if they first appear as dummy arguments 
in a FUNCTION, SUBROUTINE, or ENTRY statement. This is a valid example 

SUBROUTINE SUB (X,Y,Z,I) 



ENTRY SUB1 (A,B) 

C = A+B 

Example 2 : 

C alling Program 

CALL SUB1 (A,B,C,D,E,F) 
CALL SUB2<G,£10,£20) 



Subprogram 

SUBROUTINE SUB1 t U,V, W, X, Y, Z> 

RETURN 

ENTRY SUB2 <T,*,*) 

T=V*W*U 

ENTRY SUB3 (*,*) 

T=T+X*Y**Z 
50 IF (U-X) 100, 200, 300 
100 RETURN 1 
200 RETURN 2 
300 RETURN 

END 



CALL SUB3U10,*20) 
5 Y =A+B 



10 Y = C+D 

20 Y = E+F 



Explanation : A call to SUBl merely performs initialization. Subsequent 
calls to SUB2 and SUB3 result in execution of different sections of the 
subroutine SUBl. Then, depending upon the result of the arithmetic IF 
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at statement 50, return is made to the calling program at statement 10, 
20, or the next executable statement following the CALL. 

Additional Rules for using ENTRY 

1. A CALL may only change the value of explicit arguments (or implicit 

arguments in COMMON), It cannot affect the value of those that 
were initialized by some previous CALL. 

2« If a name is identified as a dummy argument only by its appearance 
in a given ENTRY statement, no use of that dummy argument may 
appear in statements preceding (physically) the ENTRY statement. 

3* The appearance of an ENTRY statement does not alter the rules 
regarding the placement of Statement Functions in subprograms. 

4 „ If new dimensions for an adjustable dimension array are to be 

passed to a subprogram with an ENTRY, the array name must appear in 
the argument list of the ENTRY* 

The EXTERNAL Statement 



| General Form [ 

| EXTERNAL a, b, £,...» | 

I I 

( where a,b,Cf.„ are names of subprograms that are used as argu- | 
j ments in other subprograms j 

If a FORTRAN-supplied in-line function is used in an EXTERNAL state- 
ment, it is not expanded in-line ; the function is assumed to be part of 
a library. (The FORTRAN-supplied in-line and out-of-Iine functions are 

given in Appendix D.») 

The name of any subprogram that is used as an argument in another 
subprogram must appear in an EXTERNAL statement. For example, assume 
that SUB and MULT are subprogram names in the following statements. 



hx ample 1 

Calling Program 

EXTERNAL MULT 



Subprogram 

SUBROUTINE SUB C X, ¥ # Z) 
IF CX) 4,6,6 
D = Y (X,Z**2> 



CALL SUB (A, MULT,C> 



RETURN 

END 



Explanation ; The subprogram name MULT is used as an argument in the 
subprogram SUB. The subprogram name MJLT is passed to the dummy vari- 
able Y; the variables A and C are passed to the dummy variables X and Z. 

The subprogram MULT will be called and executed only if the value of A 
is negative. 
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Example 2 



CALL SUB (A, B, MULT (C,D),37) 



Explanation : An EXTERNAL statement is not required because the subpro- 
gram named MULT is not an argument; it is executed first and the result 
Decomes the argument. 

BLOCK DATA SUBPROGRAM 

To enter data into a COMMON block, a separate subprogram must be 
written. This separate subprogram contains only the DATA, COMMON, 
DIMENSION, EQUIVALENCE, and Type statements associated with vhe data 
being defined. Data may be entered into labeled (named), but not unla- 
beled, COMMON by the BLOCK DATA subprogram. 



( ._, . , . — . _. . ._ 1 

| General Form | 

| BLOCK DATA 



END 



1. The BLOCK DATA subprogram may not contain any executable 
statements. 

2. Hie first statement of this subprogram must be the BLOCK DATA 
statement . 

3. All elements of a COMMON block must be listed in the COMMON state- 
ment, even though they do not all appear in the DATA statement. 
For example, the variable A in the COMMON statement below does not 
appear in the DATA statement: 

BLOCK DATA 

COMMON/ELN/C,A,B/RMG/Z ,Y 

REAL B(U>/1. 0,1. 2,2*1. 3/,Z*8( 3) /3* 7. 64980825D0/ 

COMPLEX C/( 2.4,3.769)/ 

END 

4. Data may be entered into more than one COMMON block in a single 
BLOCK DATA subprogram. 

5. No element may have more than one initial value assigned in the 
same program. 
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APPENDIX A: FORTRAN COMPARISON 



This appendix contains a description of the differences in the FOR- 
TRAN language supported by IBM OS and OS/VS r and by the IBM Time Sharing 
System. The FORTRAN language for IBM OS and OS/VS is described in IBM 
FORTRAN IV Language , GC28-6515. 

1 • Extensions 

TSS — Does not allow generalized subscripts and direct access I/O 

statements; no list-directed I/O; no free format input, 
OS and OS/VS — Allows the above. 

2. Call by Value 

TSS -- Treats all arguments as call-by-name whether or not they are 

enclosed in slashes. 
OS and OS/VS — Treats arguments not enclosed in slashes, and not 

declared as an array, as call -by- value. 

3 • Dummy Arguments 

TSS -- Dummy arguments may not appear in any statement until 
defined as such in an ENTRY, SUBROUTINE, or FUNCTION 
statement . 

OS and OS/VS — Restriction holds only for executable statements. 

4. ENTRY in FUNCTION Subprograms 

TSS -- The name of a FUNCTION subprogram must be used to return the 
value of the function, even though entry was made through an 
ENTRY statement. 

OS and OS/VS — The ENTRY name may be used to return the value of 
the function. 
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APPENDIX Bi SOURCE PROGRAM CHARACTERS 






Alphabetic 
Characters 



A 

B 
C 
D 
E 
F 
G 
H 
I 
J 
K 
L 
M 
N 
O 
P 

Q 
R 
S 
T 
U 
V 
W 
X 
Y 
Z 
$♦ 



EBCDIC or BCD 
Card Punches 



12-1 

12-2 

12-3 

12-4 

12-5 

12-6 

12-7 

12-8 

12-9 

11-1 

11-2 

11-3 

11-4 

11-5 

11-6 

11-7 

11-8 

11-9 

0-2 

0-3 

0-4 

0-5 

0-6 

0-7 

0-8 

0-9 

11-3-8 



Numeric 
Characters 




1 
2 
3 
4 
5 
6 
7 
8 
9 



Special 
Characters 



) 

♦ 

, (comma) 
( ♦♦ 

• (apostrophe) 
blank 

j . . 



EBCDIC or BCD 
Card Punches 


1 
2 
3 
4 
5 
6 
7 
8 
9 



EBCDIC 
Card Punches 



12-6-8 
11 
0-1 
6-8 

12-3-8 
11-5-8 
11-4-8 
0-3-8 
12-5-8 
5-8 
(no punch) 



BCDIC 
Card Punches 

4 



12 
11 
0-1 
3-8 

12-3-8 
12-1-8 
11-4-8 
0-3-8 
0-4-8 
4-8 
(no punch) 



-X 



-4 



Source programs are coded in either BCD or EBCDIC character codes; 
mixing the two, however, is not allowed. 

♦Considered an alphabetic character in EBCDIC only, 
♦♦Considered a special character in EBCDIC only. 
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APPENDIX C: OTHER FORTRAN STATEMENTS ACCEPTED BY TSS FORTRAN IV 



This appendix describes features of previously iiqp lemented FORTRAN IV 
languages that are incorporated into the IBM Time Sharing System FORTRAN 
IV language. The inclusion of these language facilities allows existing 
FORTRAN programs to be recompiled for use in IBM Time Sharing System 

with little or no reprogramming, 

READ Statement 



I General Form | 

READ b, list 

where b, is the statement number or array name of the FORMAT sta- 
tement describing the data 

list is a series of variable or array names, separated by 
commas, which may be indexed and incremented; they specify 
the number of items to be read and the storage locations 
into which the data is placed 

This statement causes data to be read from the data set associated 
with the system input, 

PUNCH Statement 



j General Form j 

PUNCH b, list 

where b is the statement number or array name of the FORMAT state- 
ment describing the data 

list is a series of variable or array names, separated by 
commas, which may be indexed and incremented; they specify 
the number of items to be written and the storage locations 
from which the data is taken 



Hie PUNCH statement causes data to be written in the data set asso- 
ciated with the system output. 

PRIN1 Statement 

| General Form j 

^„ _ . f 

PRINT b, list 

where b is the statement number or array name of the FORMAT state- 
ment describing the data 

list is a series of variable or array names, separated by 
commas, which may be indexed and incremented; they specify 
the number of items to be written and the locations in 
storage from which the data is taken 



8^ 



The PRINT statement causes data to be written in the data set asso- 
ciated with the system output, 

DATA Initialization Statement 



j General Form | 

data v ll ... P vyi l *d 1 ,.. M i n »d n / l v n + 1M ..,v /i n *i*dv + i.--#i *§, /.--. 

where y 1# ... # v are variables, subscripted variables (in which 

case, the subscripts roust be integer constants), or array 

names 

d s# ».,,d are values representing integer, real, complex, 
logical, or literal hexadecimal data constants 

ii,...,!. represent, unsigned integer constants indicating 

the number of consecutive variables that are to be assigned 
the value of ^ ,...,d 



A data initialization statement is used to define initial values of 
variables and arrays., There must be a one-for-one correspondence 
between these variables (i.e., v L , • . • ,v ) and the data constants (i.e., 
t| 1# . . . ,d ). 

E xample 1 : 

DIMENSION D(5,10) 

DATA A, B, C/5. 0,6. 1,7. 3/, D/25*l. 0,25*2.0/ 

Explanation : The DATA statement indicates that the variables A, B r and 

C are to be initialized to the values 5.0, 6,1, and 7*3, respectively. 
Also, the statement specifies that the first 25 variables in the array D 

are to be initialized to the value 1,0, and the second 25 to the value 
2.0. 

Example 2 : 

DIMENSION A(5), B(3,3), L (4) 

DATA A/ 5* 1.0/, B/9*2.0/, L/4*.TFUE./ f C/ f F0URV 

Explanation : The DATA statement specifies that all the variables in the 
arrays A and B are to be initialized to the values 1.0 and 2.0, respec- 
tively. All the logical variables in the array L are initialized to the 
value .TRUE., The letters T and F may be used as an abbreviation for 

.TRUE. and .FALSE. , respectively. Also, the variable C is initialized 

with the literal data constant FOUR. 

An initially defined variable, or variable of an array, may not be in 
olank common; however, in a labeled common block, they may be initially 
defined only in a block data subprogram (see m SUB PROG RAMS" ) . 

DOUBLE PRECISION Statement 



| General Form 1 

, _„„__„__— .__. ..__ „__„ H 

j DOUBLE PRECISION a,b,c # ... ( 

I I 

| where f*,b,c,... are variable names that may be dimensioned in the | 

j statement, or function names j 
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The DOUBLE PRECISION statement explicitly specifies that the 
variables a f b,c f ... are of type double precision. This statement over- 
rides any specification of a variable made by either the predefined con- 
vention or the IMPLICIT statement. This specification is identical to 
that of type REAL*8. 

Also, FUNCTION subprograms may be typed double precision, in this way 

DOUBLE PRECISION FUNCTION name (a , ,a 2 f a 3 , . . . f a n ) 
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APPENDIX D: FORTRAN SUPPLIED SUBPROGRAMS 



The FORTRAN supplied subprograms are of either of two types: mathe- 
matical subprograms and service subprograms. The mathematical subpro- 
grams correspond to a FUNCTION subprogram; the service subprograms 
correspond to a SUBROUTINE subprogram- Appendix D lists the in-line and 
out-of-line mathematical FUNCTION subprograms. An in-line subprogram is 
inserted by the FORTRAN compiler at any point in the program where the 
function is referenced* An out-of-line subprogram is located on a 
library. A detailed description of out-of-line mathematical subprograms 
and service subprograms is given in FORTRAN IV Library Subprograms . 



MATHEMATICAL SUBPROGRAMS 

All functions are used as described in the section "FUNCTION Subpro- 
grams m — i.e., A = AMOD(X i# X 2 ) # where A is the value and X t and X 2 are 
the arguments. 

Table 4. Mathematical function subprograms Cpart 1 of 3) 



Common Logarithm 



Arcsine 

Arccosine 

| 

Arctangent 



Function 
Exponential 



Natural Logarithm 



— I- 



ARJIN 
DARSIN 

+ 

ARCOS 
D ABC OS 



- + 



Trigonometr ic 
Sine 

(Argument in 
radians) 

r 

Trigonometric 
Cosine 
(Argument in 
radians) 



Trigonometric 
Tangent 






Trigonometric 
Cotangent 



Square Root 



Hyperbolic 



Name 

EXP 
DEXP 
CEXP 
CDEXP 



ALOG 
DLOG 
CLOG 
CDLOG 



ALOG 10 
DLOG10 



Def i nit i on/Us age 

earg 
earg 
earg 
earg 
A=EXP(Xi> 

In (Arg) 
in (Arg) 
In (Arg) 
In (Arg) 
A=ALOG(X,) 

log l0 (Arg) 
log l0 (Arg) 
A^ALOGKMXi) 



AT AN 
ATAN2 

DATAN 
DAT AN 2 



1 J£ 



SIN 
DSIN 
CSIN 
CDSIN 



COS 
DCOS 
CCOS 
CDCOS 



TAN 
DTAN 



COTAN 
DCOTAN 



SQRT 
DSQRT 
CSQRT 
CDSQJRT 



SINH 
DSINH 



arcs in (Arg) 
A=ARSIN(X t ) 

arcos (Arg) 
A=ARCOS(X,) 



4 



4 — 



arctan (Arg) 
arc tan CArgj/Arga) 
arctan (Arg) 
arctan ihrq l /hrq 2 } 
A=ATAN(X, ) 



sin (Arg) 
sin (Arg) 
sin(Arg) 
sin (Arg) 
A=SIN(X,) 



cos (Arg) 
cos (Arg) 
cos (Arg) 
cos (Arg) 
A=COS(X!) 

tan (Arg) 
A=«TAN(X, ) 



co tan (Arg) 
A=COTAN(X t ) 



(Arg) 
(Arg) 
( Ar a ) 

A=SQKT(X! ) 



sinh (Arg) 
A=SINH(Xj ) 



In-Line (I) 
Out-of-Line (O) 



-r 

1 
1 

4 

1 



o 

o 

O 
O 



4 4 l_ 



o 

o 
o 
o 



o 
o 
o 
o 



o 
o 

o 
o 



No. Of 
Arg. 



Type of [ 
Arguments [Function 
+ 

Real *«t (Real ♦<* 

Real *8 [Real *8 

Complex *8 [Complex *8 
Complex *16 (Complex *16 
I 

Real *<* [Real *M 

Real *8 [Real *8 

Complex *8 [Complex *8 

Complex *16 [Complex »16 

[ 

4 _ 

Real •<* [Real *<* 

Real »8 [Real *8 



Real *4 
Real *8 



Real *H 

Real *8 

Real **» 

Real *<* 

Real *8 

Real *8 



[Real *4 

[Real »8 

H 

[Real *U 
| Real *8 



[Real *« 
[Real ♦ <* 
(Real *8 
| Real *8 



Real *4 
Real »8 
Complex 
Complex 



_4_- 



[Real *4 
[Real *8 
♦8 (Complex *8 
♦16 | Complex *16 
I 

Real *** (Real *4 

Real *8 (Real *8 

Complex *8 (Complex *8 

Complex *16 [Complex *16 



Real *4 
Real *8 



Real ♦<* 
Real *8 



H 

[Real *4 
|Real *8 

-4 

[Real *4 
(Real *8 



Real ♦<* [Real *** 

Real *8 [Real *8 

Complex *8 [Complex *8 

Complex *16 | Complex *16 



Real ♦*» 
Real *8 



H * 

[Real •** 
|Real *8 
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Table 4. Mathematical function subprograms Cpart 2 of 3) 



Complemented 

Error Function 

Gamma 



Log -gamma 



Remainder ing 



Absolute value 



Hyperbolic 
Cosine 



-4 



Hyperbolic 
Tangent 



Error Function 



Name 

f— - 

COSH 
DCOSH 

+ — 

TANH 
DTANH 



ERF 

DERF 



ER FC 
DERFC 

4- 

GAMMA 
DGAMMA 



ALGAMA 
DLGAMA 

4 

MOD 

AMOD 

DMOD 



IABS 
ABS 
DABS 
| 

CABS 
CDABS 



Truncation 



Largest value 



Smallest value 



Transfer of sign 



. + , 



Positive 
difference 

Obtaining most 
significant part 

of a Heal *8 
argument 

Obtain real 
part of complex 
arg oment 



A=ERF(X { > 

1-erf <x) 
A=ERFC(x,> 



INT 



AINT 
ID INT 

AMAXG 

AMAX1 

MAXO 

MAX1 

DMAX1 

AMINO 
AM INI 

MI NO 
fci INI 
DMIN1 

FLOAT 
D FLOAT 



IFIX 
HFIX 



ISIGN 
DSIGN 



DIM 
IDIM 



REAL 



Definition/Usage 



cosh (Arg) 
A=COSH(X l > 

tanh(Arg) 
tanh (Arg) 
A=TANH(X, ) 



In-Line (I) 
Out-of-Line CO) 

I O 



/ 



~u 2 

e du 



x-1 -u 
u e du 



ANS=GAMMA(X, > 



log c r(x) 
A=ALGAMA(x, ) 

Arg , (mod Arg ? ) 
A=MOD(X , ,X ? ) 



| Arg | 
A=ABS(Xi> 



jVx ' +x 2 f for Xj+x, 
A=CABS(X.) 



Sign of Arg times 
largest integer 
<|ArgI 

I=INT(X, ) 

Max (Argj # Arg , r . . . ) 



A=AMAXO(Xj ,X 2 ,.. X n ) 
Min (Arg, § Arg ? ,. .. ) 

A=AMINO(X, ,X ?# ..X n ) 

Convert from 
integer to real 
A=FLOAT(X, ) 



Convert from 
real to integer 
I^FIXCX^ 

Sign of Argj times 
1fcrg,| 

A=SIGN(X l ,X 2 ) 

Arg t -Min (Arg , *Arg 2 ) 
A=DIM(X, ,X ? ) 



A=SNGL(X,) 



A=SNGLIX, ) 



No. of 

Arg. 



4 

1 

1 



4 



1 

1 



— — V 



>2 
>2 

>2 

>2 

>2 
22 

2.2 

>2 
22 

1 
1 



1 



-— 4 
i 



Type of 
Arguments 

Real *4 

Real *8 

Real *a 
Real *8 



Real *a 

Real *8 



Real *tJ 
Real *8 



Real *4 

Real *8 



Real *4 

Real *8 



Int eg er * H 
Real *U 
Real *8 

Integer *4 
Real *U 

Real *8 



Complex *8 
Complex *16 



Real *4 



Real *H 
Real *8 

Integer *H 
Real *4 
Integer *1 
Real *U 
Real *8 

Integer *4 
Real *** 
Integer *4 
Real *a 
Real *8 

Integer *** 
Integer *<* 



Real *U 
Real *ti 



Real ♦«* 



Integer **4 
Real *8 



Real *4 

Integer *H 



Complex *8 



Function 



Real **4 
Real *8 



Real *a 
Real *8 



Real *U 
Real *8 



Real *U 
Real *8 



Real *U 
Real *8 



Real **4 
Real *8 



Integer *4 
Real *a 

Real *8 

Integer *U 

Real *U 
Real *8 



Real *4 
Real *8 



Integer *4 



Real *M 

Integer **+ 

Real *4 
Real ♦«* 

Integer *4 
Integer *4 
Real *8 

Real *U 
Real *U 

Integer *** 
Integer *ft 
Real *8 

Real *4 
Real *8 



Integer *4 

Integer *2 



Integer *4 
Real *8 

Real *4 

Integer *U 
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Table 4. Mathematical function subprograms (part 3 of 3) 



T . — . T __. . r — _. — . — 

No. of J Type of | 

Arg. j Arguments j Function 

1 {Complex *8 |Real »U 



I 

| Function 

J Obtain imaginary 

} part of complex 
[argument 

j Express a Real 
|**4 argument in 
(Real *8 form 

[Express two real 
| arguments in oom~ 
Iplex form 

j Obtain conjugate 
| of a complex 
|argument 



|- 



Name (Definition 



AIM AG | 
I 



In-Line (I) 
| Out- of- Line CO) 

I 



|A=AIMAGCX, ) 



DBLE 






1 | Real *a 



|A=DBLE(3C, ) 

CMPLX lC=Arg, +ihrg 2 
DCMPLX | 

|A=CMPLXCX ,X 2 ) 



I 

I 



CONJG |C=X-iY 
DCONJGfFor Arg==X+iY 
|A=CONJG£Xi ) 






2 JReal *U |Complex *8 

2 fReai *8 | Complex *16 

I I 

1 [Complex *8 [Complex *8 

1 [Complex *tb | Complex *16 

I I 



SERVICE SUBPROGRAMS 

MACHINE INDICATOR TESTS SUBPROGRAMS 

In the list of pseudo machine -indicator test subroutines below, 
assume that i is an integer expression and that j is an integer vari- 
able. These subroutines are referred to by CALL statements. 

SLITE (i) : If i = 0, all sense lights will be turned off. If i = 1, 2 f 
3, or 4, the corresponding sense light will be turned on, 

SLITET (1,1) : Sense light i (equal to l r 2, 3 r or 4) will be tested and 
turned off. The variable j will be set to 1 if i was on, or j will be 
set to 2 if i was off. 

Example : Assume that the program is to continue if sense light i is on 
and the results are to be written if sense light i is off. This can be 
done by using the logical IF statement or a computed GO TO statement: 



CALL SLITET (3, KEN) 
GO TO C6, 17) ,KEN 

17 WRITE (3, 26) CANS CK) 

6 CONTINUE 



K=l, 10) 



Explanation ; When the statement CALL SLITETC 3, KEN) is executed, the 
variable KEN is assigned the value 1 or 2 depending on whether sense 
light 3 is on or off, (and the sense light is turned off). If KEN is 1, 
statement 6 is executed next? if KEN is 2, statement 17 is executed. 

QVERFL Cj) : j is set to 1 if a floating-point overflow condition 
exists, i.e., if the result of an arithmetic operation is greater than 
16 63 ; j is set to 2 if neither an overflow condition or underflow condi- 
tion exists; j is set to 3 if floating-point underflow condition exists, 
i.e., if the result of an arithmetic operation is less than 16-° 3 . The 
machine is left in a no-overflow condition. If a sequence of operations 
caused both overflow and underflow to occur, the value of j returned 
represents whichever of these two conditions occurred last. 

DVCHK Cj) : If the divide c eck indicator is on, j is set to 1 and the 
divide check indicator is turned off; if the divide check indicator is 
off, j is set to 2. 



Appendix D: Fortran Supplied Su bprograms 89 



THE EXIT, DUMP, AND PPUMP SUBPROGRAMS 
EXIT Subprogram 

A CALL to the EXIT subprogram terminates the execution of the object 

program. 

PUMP Subprogram 

A CALL to the DUMP subprogram 

CALL DUMP (A 1 ,B I ,F lf ...,A n ,B n ,F„> 

causes the indicated limits of storage to be duinped and execution to be 
terminated. 

1. A and B are variable data names that indicate the limits of storage 
to be dumped; either A or B may represent upper or lower limits. 

2. F n is an integer indicating the dump format desired 



n = 


hexadecimal 


1 

2 
3 
4 
5 


logical *1 
logical *1 
integer *2 
integer *** 
real *4 


6 


real *8 


7 
8 
9 


complex * 8 
complex*16 
literal character 



3. If the argument F n is omitted, it is assumed to be equal to 0, and 
the dump will be hexadecimal. 

4. The arguments A and B should be in the same program (main program 
or subprogram) or same common block. 

PPUMP Subprogram 

A CALL to the PPUMP subprogram 

CALL PPUMP (A, f B, ,F l ,.. . ,A n ,B n ,F n ) 

causes the indicated limits of storage to be dumped and control to be 
returned to the calling program. 



**Q 



APPENDIX E; EXAMPLES OF FORTRAN- WRITTEN PROGRAMS 



EXAMPLE PROGRAM 1 

Example program 1 (Figure 2) is designed to find all the prime num- 
bers between 1 and 1000, A prime number is an integer that cannot be 
evenly divided by any integer except itself and 1, Thus l r 2, 3, 5, 1, 
11 # ... are prime numbers. The number 9 f for example, is not a prime 
number, since it can be divided evenly by 3. 



IBM 



SAMPLE PROGRAM I 



6/66 



8! 
ji 
101 

3 
107 
103i 
104 ! 
105 r 
106' 

1' 
107 

5" 

2 + 

108; 

'I 

109 



PRIMiE NUMBER PpOBLEJM ] | 

WRITE (6>B) 1 1 

FORWJAT (5i2H FOLLOWING I S| A LIST OF 
.:i?XviHl/ij9XilH2/19X*lH3)! 

!l«5 I 1 . ' • ' 

Mi!; : 

jA-SQRT(A): ! 

<J«A 1 [ 

!D0 1' K = 3iJi2 ' 

!l • i /k * , 

IF(L'*K- I )'i i2>MJ 

'CONTINUE ' ! 

Write (6vs)r ! 

FORMAT ( 120) ] 
'I » I +2 ! 

'iFdjwtf-r 1 ) 7 -> 4V3 

Write (&i9) 

iformat (1«h program errojr) 

Write (6->%) 

FORMAT (3 IN THIS IS 
STOP 

END 



prijaJe NUM&ERS 






IF ROM ti TO 



THE ]END OJF THE 



PROGRAM) 



■J~± 



1000/ 



Figure 2. Example Program 1 



EXAMPLE PROGRAM 2 



The n points Cx , y ) are to be used to fit an m degree polynomial by 

the least -squares method. 



y = a + ajX + a 2 x a ♦ 



4- a xm 



To obtain the coefficients a # a+,... # a 9 it is necessary to solve 
these normal equations 

CI) W a + M l ck l + ... ♦ Wmam = Z 

C2) W t a + W 2 a, ♦ ... ♦ Wm+ ^ = Z x 



(m+1) 



Wma + Wm+ j a , + ... * W 2 raam = Zm 



where 



W = n 



Z = 2 yi 
i=l 
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Wi » 2 xi 



Z. = 2 yixi 



n 
W 2 = 2 xi* 
i=l 



n 



Z 2 = 2 yixi- 
i=l 



n 
Zm * 2 yixim 
i=l 



W 2 m = 2 xi a m 
i=l 



After the Ws and Zs have been computed, the normal equations are 
solved by the method of elimination, which is illustrated by the follow- 
ing solution of the normal equations for a second degree polynomial (m = 

2) . 

(1) W a ♦ W^! + W 2 a 2 = Z 

(2> W^o ♦ W 2 aj ♦ W 3 a 2 = Z x 

(3) W 2 a ♦ W 3 a x ♦ W a 2 = Z 2 
Hie forward solution is 

1. Divide equation (1) by W . 

2. Multiply the equation resulting from step 1 by W,, and subtract 
from equation ( 2) . 

3. Multiply the equation resulting from step 1 by W 2 , and subtract 
from equation (3), 

The resulting equations are 

(4) a ♦ b l2 a, + b l3 a 2 = b 14 

(5) b 22 a, ♦ b 23 a 2 = b 24 

(6) b 32 ai ♦ b 33 a 2 — b 34 
where 

b l2 = W v /W 0# b 13 = W 2 /W , b, = Z /W 

b 22 = W 2 -b l2 Wi , b 23 = W 3 -b u Wi , b 24 = Z^bj Mm 

b ]z = W 3 -b 12 W 2 , b 33 =W 4 -b 13 W 2 4, b 34 = Z 2 ~bj W 24 

Steps 1 and 2 are repeated, using equations (5) and (6), with b 8a and 
fa 32 instead of W and W x . The resulting equations are 



(7) 
(8) 



l i 



^ i * *■* 2 3 ** 2 2= C 2 4 
C 33 a 2 = C 34 



92 



where: 



= b ,/b 



, c 24 = b 2 /b 2 



c 33 - b 33 -c 23 b 32 , c 34 = b 3 -c 24 b. 

The backward solution is 

(9) a 2 = c 34 /c 33 

(10) a j = c 24 -c 23 a 2 

(11) a = b 14 -b l2 a l -b 13 a 2 



from equation (8) 
from equation (7) 
from equation (**) 



Figure 3 is a possible FORTRAN program for carrying out the calcula- 
tions for the case: n = 100, m < 10. W 0f W lf W 2 , . .., W 2 m are stored 

in W(l), W(2), W(3) f . .,, W(2M+1), respectively, Z , Z lf Z 2 , , Zm 

are stored in Z (1), Z (2), Z (3), •••, Z(M+1) # respectively. 



IBM 



MTRAN CsrfiRf r.rrr 



5AMPLE PROGRAM 1 



p/66' 



i : (Jvlrff 



jREALj X(l< 



(FORM 



AT (I 



FORMAT (5 



I 
13| 



READ 

|LW - 

LB = 
LZ = 

bo 5 
wt( J) 
w'l) 
bo 6 
z(J> 
bo i 
jp = 

|2(i) 
DO I 
p = 

W( J) 
Z(J) 



p = 



(5il 
2*M+ 
M«-2 
M+t 

J*2 

= 0. 

* N 

J*i 

« 0. 

|6 I- 



2( 

3 J = 

X(I) 

- W( 

- Z( 



)>y(|i00)iw(2i)|->z(ii 

2»I3/IUFH*..7) } 
El 5.61) 

) min!i(x(i|) >yu 

i I 



DO 16 J 



X(I)*P 



iLW 





iLZ 




t)+Y(D 
2»L2 I 



HP 



J)+P ! 
J)+Y(;I)*P 



) i I « i 



)>A<1 



>N) 



i)»B(iiiri2[) 



Figure 3. Example Program 2 (part 1 of 3) 
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ISM rOHTRAK Cattnf T,rm ■ i" '. 


SAMPLE PROGRAM 2 




















... 2^T 




_. _ 


1 ■ 6/66 
























'0»'l*N i'A'fMtMl 


•"rrj'; " 


16 
! 17 








W( J) 
DO 2 
DO 2 


- W(|J)+P 


1 




















■- 




F = 

K+I 

I) - 
2 K* 


1*12 
1>L2 

WCJ-1 
1>L2 


) 














H 










1 M = 

1 2« BCKi 

I ! !DO 2 

22 B(Ki 

2 3J DO 3 

1 ! JDIVB 

! ,DO 2 

26 B(Li 

: ! IF ( 

28 DO 3 









































— 




LB) - 2(K) 
1 L»i^L2 

= B(LiL) 
6 J-liLB, 
J) - |B(LiJ 

L + i 
Il-LB) 28? 
i I = 1 1 » L 2 
TB = B(IiL 


J 

















J 










J -■--( 


}/_DIV 
33*33 

). _ 


B 

























.._. 






















I 


iFMUL 
















_ 


h l J4,LB 


)-B(L 
BT""" 


>J)*F 


MULTB 



















1 3i{ |B(I»J) * |B( hJ 














33 A(LZ) = B 


(LZ^L 


















' i 1 • F z 1 




. _-_ _. 













(• ■• 




3 5, SIGMjA = 0-P 














L ... . I 


f P° S 


7 J» 


IiLZ 















4 
















_L L-- 




L _ _ 










I 


i 




"~1 














■• ■. •. . - . .. «i 





Figure 3. Example Program 2 (part 2 of 3) 
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^J 
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FO»T«AN SIMEmCNI 




£ 


37 
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B f I " 1 


>J >*A 


(J) 


• ' [ 




..,,.. 














i- = i-i r 








1 














1 




AfT) 1 


- B(|I»LB) 


-ugh 

>35 

(IM 


/L ... 
= 1>Z) 


— 


J 












..__ _ , 


> ■■ 


4i 


5f ( 

4 WRIT 
STOP 










■ ■ ( - - j 




i- 




END 




























































I 


... 


















^ 


_ 


_ _ . 




___ 








1" 














































































i i . 












- - 








•• •— — " 


i 




| 




























t 


.i_ 































!— - 

-I — 





























I ■ i 
- 

f--4 


-~ 

















L - 


i 



Figure 3. Example Program 2 (part 3 of 3) 
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The elements of the W array, except W(l), are set equal to zero? W(l) 
is set equal to N. For each value of I , X and Y are selected. The 
powers of X are computed and accumulated in the correct W counters. 
The powers of X are multiplied by Y , and the products are accumulated 
in the correct Z counters. To save machine time when the object program 
is being run, the previously computed power of X is used when computing 
the next power of X . Note the use of variables as index parameters. 
By the time control has passed to statement 17, the counters have been 
set to 



W(l) = N 



Z(l> = 



N 

2 Y s 
1=1 



W(2) = 2 Xj 
1=1 



Z(2> = 



N 

2 YiXj 
1=1 



W(3) = 2 X a 
1=1 



Z(3> = 



N 

2 YjXj* 
1=1 



Z(M+1) =2 y. XjM 
1=1 



W(2M*1) = 2 X 

i 

1=1 



2 M 



By the time control has passed to statement 23 r the values of W , Vl { , 
. .., W 2 m+ L have been placed in the storage locations corresponding to 
columns 1 through M + 1, rows 1 through M + 1, of the B array, and the 
values of Z , Z l9 . .., Zm have been stored in the locations correspond- 
ing to the column of the B array. For example, for the illustrative 
problem (M = 2), columns 1 through 4, rows 1 through 3, of the B array 
would be set to these computed values 



W n 



w, 



w s 






w 2 
w, 



W 4 



z 2 



This matrix represents equations CI), (2), and (3), the normal equa- 
tions for M = 2. 



The forward solution, which results in equ 
the illustrative problem, is carried out by s 
By the time control has passed to statement 3 
Al terms in the M ♦ 1 equations, which would 
1 at ions, have replaced the contents of the lo 
columns 1 through M ♦ 1, rows 1 through M + 1 
constants on the right-hand side of the equat 
tents of the locations corresponding to colum 
11 ♦ 1, of the B array. For the illustrative 
4, rows 1 through 3, of the B array would be 
values 



ations <U) , (7), and (8) in 
tatements 23 through 31 . 
3, the coefficients of the 
be obtained in hand calcu- 
cations corresponding to 
, of the B array, and the 
ions have replaced the con- 
n M ♦ 2, rows 1 through 
problem, columns 1 through 
set to these computed 
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1 


b, 2 


b n 





1 


C 23 








c 33 



b x 



This matrix represents equations (4) , (7) , and (8). 

The backward solution, which results in equations (9), (10), and (11) 
in the illustrative problem, is carried out by statements 33 through 40 , 
By the time control has passed to statement 41, which prints the values 
of the A9 terms, the values of the (M + 1) *Aj terms have been stored in 
the M + 1 locations for the A array. For the illustrative problem, the 
A array would contain these computed values for a 2r a x f and a 

Location Contents 

A(3) c 34 /c 33 

A(2) c 2 4 -c 23 a 2 

A(l) b l4 -b 12 aj-b l3 a 2 

The resulting values of the AI terms are then printed according to 
the FORMAT specification in statement 2. 



96 



INDEX 



Where more than reference is given, the 
major reference is first. 



A format code 49-51 

AB5 88 

absolute dimensions of an array 60-61 

absolute value 88 

addition 14 

adjustable dimensions of an array 60-62,1 

AIMAG 89 

A INT 8 8 

ALGAMA 88 

ALOG 8 7 

ALOG10 87 

alphabetic characters 83 

alphameric data format code 49-51 

AMAXO 88 



88 
88 
88 
88 
18 



87 



82 
52,83 



AMAX1 

AMINO 

AMI Ml 

AMOD 

• AND. 

arccosi ne 

ARCOS 87 

a re sine 87 

arctangent 87 

arguments, dummy 

literal constant 

arithmetic assignment statement 

arithmetic condition 17 

arithmetic expressions 13-17 

arithmetic IK statement 24-25 

arithmetic operators 14,15 

arrangement of arrays in storage 

arrays 10-13 

absolute dimensions of 60-61 
adjustable dimensions of 60-62,1 
name of an array 33,66-68 
size 12 
subscripts 10-11 

ARSIN 87 

ASSIGN statement 23-24 

assigned GO TO statement 23-24 

ATAN 8 7 

ATAN2 87 

attribute control, variable 1 



BACKSPACE 56,31 
BCDIC card punches 83 
beginning of a data group 33 
beginni ng-of-data transfer code 53 
blank common area 64-65 
blank output records 42 
BLOCK DATA subprogram 81 
boundary alignment 65-66 



20-21,2 



12-13 



card input 2 


carriage control 


ccos 


87 


CDABS 


88 


CDCOS 


87 


CDEXP 


87 


CDLOG 


87 


CDSIN 


87 


CDSQRT 87 


CEXP 


87 



55-56 



CABS 8 8 

CALL statement 76-77 

call-by -value 82 



characters in program 83 

CLOG 87 

CMPLX 89 

coding FORTRAN statements 2 

comments 2 

COMMON 6 2-66,57,12 

blank 64-65 

implicit arguments 8 2 

named common area 64-65 

statement 6 2 
common logarithm 87 
compiler 1 

complemented error function 8 8 
complex 

argument 

obtaining conjugate 89 
obtaining imaginary part 89 
obtaining real part 88 

COMPLEX 57-59,15 

constant 5 

real arguments expressed in complex 
form 8 9 

statement 57-59 
CONJC 89 
constants 3-8 
continuation line 2 
CONTINUE statement 2 9 
control statements 22-3 0,2 
COS 87 
COSH 8 8 
COTAN 87 
CSIN 87 
CSQRT 87 



D format codes 48 f 46 

DABS 8 8 

D ARCOS 87 

DARSIN 87 

data initialization 85,6,59 

data set 31 

DATAN 87 

DAT AN 2 87 

DBLE 8 9 

DCMPLX 89 

DCONJC 89 

DCOS 87 

DCOSH 88 

DCOTAN 87 

decimal places 4 6,48 

decimal point 48 

DERF 8 8 



Index 97 



positive 88 



DERFC 88 
DEXP 8 7 
D FLOAT 88 
DSIN 8 7 

DGAMMA 8 8 

difference, 
DIM 88 
dimensi on 

adjustable 94 

DIMENSION 6 0- 61 , 1 2 , 57 
displaceraent r array 68 
division 14,17 
DLGAMA 88 
DLOG 8 7 
DLOG10 87 
DMA XI 88 
DMIN1 88 
DMOD 8 8 

DO loop programming 
DO statement 26-27 
DOUBLE PRECISION 85 -86 
DSIGN 88 
DSINH 87 
DSQRT 87 
DTAN 8 7 
DTANH 88 
dummy argument 
dummy array 61 
dummy variables 
DUMP subprogram 
DVCHK 89 



function 6 9 

function subprograms 
return of value 82 



72-75, 69,87 



considerations 28- 29 



71,82 

76,68 
90 



E format codes 48,46 

EBCDIC card punches 83 

elements of the language 2 

&END 33-34 

END 30,31,72 

END FILE 56, 31 

end of a data group 34 

r,NTRY statement 78-80,82 

EQUIVALENCE 66-68,57 

ERF 88 

ERFC 8 8 

ERR 31 

error function 8 8 

EXIT Subprogram 90 

EXP 87 

explicit specification 59-60,9-10,57 

exponential 14-15,87 

expressions 13-19 

EXTERNAL statement 8 0-81 



F format code 47,4 6 

fix 88 

FLOAT 88 

FORMAT 4 0-56,1, 31 

alternate Hollerith 6 

code s 4 2-53 

statement rules 40-42 
FOR IRAN 

codi ng form 3 

differences between OS 
FORTRAN 8 2 

records 41-42 

special TSS FORTRAN IV 

supplied subprograms 



or OS/VS and TSS 



features 
87-90,69 



G format code 42-46 

GAMMA 88 

general format code 42-46 

GO TO statements 22-24 
assigned GO TO 23-24 



H format code 52 
hexadecimal constants 6-7 
HFIX 8 8 

hyperbolic cosine 8 8 
hyperbolic sine 87 
hyperbolic tangent 88 



I format code 47,4 6 

IABS 8 8 

IDIM 88 

IDINT 88 

IF statement 2 4-26 

IFIX 88 

imaginary part of a complex argument 89 

IMPLICIT 57-58 

IMPLICIT specification 9 

in-line subprogram 87 

indexing I/O lists 36-37 

initialization of data 85 f 6, 59 

input data 33-34 

input/output statements 31-56 

insert blanks 52-53 

INT 88 

INTEGER 57-60,12 

integer constants 3-4 

integer data format code 4 7 f 46 

integer division 17 

integer mode 14 

I/O lists 31-37,38-39 

ISIGN 88 



keyboard input 2 



L format code 4 9 

labeled common area 64-65 

largest value 8 8 

length 

specification for variables 

total field length 44-4 5 
list 31-37,38-39 

(see also NAMEL1ST) 
literal constants 6, 1 
literal data 51-52 
literal format code 1 
log- gamma 8 8 
logical 

assignment statement 20-21 

constants 6 

expressions 17-19,13 

format code 4 9 

IF statement 25-26 

operators 18 



7 r 44-45 
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machine indicator tests subprograms 89 
mathematical subprograms 87 
MAXO 8 8 

MAX1 8 8 

M IN 8 8 

MINI 8 8 

mixed mode 1 

MOD 88 

mode 13,1 

modulus 88 

multiline lifting 43-45 

multiple entry into a FUNCTION 

subprogram 7 5 
multiplication 14 



named common area 64-65 

NAMELIST 3 2-34,38-39 
I/O 1 

name 3 2, 3 3 

natural logarithm 87 

nested DOs 28 

.NOT. 18 

numeric characters 8 3 

numeric format codes 47-49,46 

numerical constants 3-5 



conjugate of a complex 
imaginary part of complex 
real part of complex 



obtaining the 

argument 89 

obtaining the 

argument 89 

obtaining the 

argument 88 

operators 

logical 18 
relation of 17 
optional length specification of 

variables 8-9 
. OR . 18 
order of compution 19 

in arithmetic expressions 16 
in logical expressions 19 
OS/TSS FORTRAN differences 82 
out-of -line subprogram 87 
OVERFL 89 



P scale factor 54-55 

padding character 47-49 

parentheses 

in arithmetic expressions 16-17 
in logical expressions 19 

PAUSE statement 30 

PDUMP subprogram 90 

positive difference 88 

PRINT statement 84-85 

printing 55-56 

programs, sample FORTRAN 91-96 

PUNCH 84 



READ 31-37 

READ lists 31-37 

reading FORMAT statements 37 

real 

complex form, two real arguments 
expressed in 89 



constants 4-5 

mode 14 

REAL 9,57-60,88 
relational operators 17 
remaindering 8 8 

RETURN statement 88-89,77,72-73,75-76 
REWIND 56,31 

sample FORTRAN programs 91-9 6 

scale factor 54-55 

service subprogram 89, 87 

SIGN 88 

SIN 87 

SINH 87 

skip characters 52-53 

SLITE 89 

SLITET 89 

smallest value 88 

SNGL 8 8 

source program characters 83 

spacing format code 53,1 

SQRT 87 

square root 87 

standard length specification of 

variables 7 
statement 2 

control 2 2-30 

functions 70-71,69 

input/output 31-56 

numbers 2 

specification 57-68 

subprogram 69-81 
STOP statement 3 
storage specification 60-6 6 
subprogram 69-81, 87-89 

name 6 9 

service 89,87 

statements 2 
subroutine subprogram 75-81,87 
subscript 10-11, 82 
subtract! on 14 



T format code 53 

TAN 87 

TANH 88 

transfer of sign 88 

trigonometric cosine 87 

trigonometric cotangent 87 

trigonometric Sine 87 

trigonometric tangent 87 

truncation 88 

TSS FORTRAN IV 2 

special features 2 

TSS/OS, OS/VS FORTOAN differences 8 2 
type specification of the FUNCTION 

subprogram 73-74 
type statements 57-6 



value, call-by 82 
variable 8-10 

attribute control 1 

length specification 7,89 

names 8,31 

predefined specification of variable 
type 7 

types 8-9 



99 



WRITE 37-3 9,31 
WRITE lists 38-39 
writing blank lines 



X format code 52 



39 



Z format code **8-49,«*6 
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